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FOREWORD 


This is the third volume in a final report series describing 
work completed under Project NAS 9-14467. The overall objectives of 
the project are concerned with applying canopy reflectance modeling 
to signature extension tasks for wheat identification. This effort 
supports the LACIE program. Mr. T. Barnett is overall technical 
monitor of the project. Mr. M. McEwen is technical monitor for the 
Field Measurements Program. 

Volume I presents the multiplicative and additive coefficient 
matrices for a linear sun-angle correction approach. These coeffi- 
cient tables are calculated using either measured empirical canopy 
reflectance functions or model derived data. These values are then 
incorporated into an atmospheric radiation transfer model. The 
dependence of the coefficient matrices on crop stage, crop type, and 
canopy directional reflectance variations is reviewed. Finally, a 
method for inferring leaf area index, an intrinsic scene characteristic, 
from canopy reflectance is discussed. 

Volume II presents the basic data and computer programs used in 
the study, A brief review of the radiometric and geometric data 
collection procedures is also given. In particular, two recent methods 
developed by the investigators for deteraining plant geometry are 
discussed. These include the Fourier diffraction and multiple view 
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angle approach. The data compilation consists of canopy reflectance 
Leaf-Area- Indices, and leaf slope distributions for four wheat crop 
development stages at Garden City, Kansas. 



ABSTRACT OF THESIS 

EXTRACTING SCENE FEATURE VECTORS THROUGH MODELING 

The science of Remote Sensing has traditionally been involved 
with the identification of objects and materials, through varia- 
tions in electromagnetic fields. More recent developments have 
extended this research into inferrences about target status or con- 
dition. The development of relationships used in estimating scene 
status has been primarily based on empirical study. This report 
represents a different approach, in that it utilizes computer mod- 
eling in the development of a data base expressing the relationships 
under investigation. 

The specific area of research is the remote estimation of leaf 
area index of winter wheat at Finney County, Kansas. The procedure 
developed consists of three activities: 1) field measurements; 

2) model simulations; and 3) response classifications. The first 
activity is designed to identify model input parameters and develop 
a model evaluation data set. A stochastic plant canopy reflectance 
model is employed to simulate reflectance in the LANDSAT bands as a 
function of leaf area index for two phenol ogi cal stages. An atmos- 
pheric model is used to translate these surface reflectances into 
simulated satellite radiance. A divergence classifier determines 
the relative similarity between model derived spectral responses and 
those of areas with unknown leaf area index. The unknown areas are 
assigned the index associated with the closest model response. 



This research demonstrated that the SRVC canopy reflectance 
model is appropriate for wheat scenes and that broad categories 
of leaf area index can be inferred from the procedure developed in 
this study. The evaluation data set was insufficient for testing 
the procedure's accuracy in predicting specific indices. 

Other significant contributions of the study include the devel- 
opment and refinement of two field techniques for assessing leaf 
angle distribution, and the presentation of an empirical data set 
containing environm, .r..al factors, intrinsic scene parameters and 
canopy reflectance for a single area throughout two phonological 
stages. 
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1.0 INTRODUCTION 


An enduring characteristic of man is his earnest efforts to 
domesticate his environment. This process of modifying wildland 
conditions in order to benefit mankind is predicated on his ability 
to survey the earth's surface; in order to learn about its resour- 
ces and thus manage them better. Questions as to how this informa- 
tion can be derived and the mechanisms necessary to extract it have 
fostered the establishment of the field of remote sensing. 

Remote sensing is the science of acquiring information about 
material objects from a distance. This process involves two prin- 
cipal activities: 1) the procurement of physical measurements; and 

2) the translation of these data into useful information about ob- 
jects. In remote sensing, physical information may be transmitted 
to the observer either through force fields or electromagnetic 
fields; in particular, through the spectral, spatial, and temporal 
variations of these fields. Therefore, in order to derive useful 
information from these field variations, one must be able to 

- measure the variations and 

- relate these measurements to those of known 
objects or materials. 

Of the two types of fields, electromagnetic fields provide the great- 
est contemporary use. In particular, the visible and near infrared 
portions of the electromagnetic spectrum dominate the science of 
remote sensing. The remainder of this paper is confined to fields 
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of this type, and procedures for translating these data into infor- 
mation about the status of the target. 

Previous development of the methodologies for relating physical 
measurements to information, have been principally based on empiri- 
cal studies. These investigations correlate measured changes in 
one variable, or set of variables, to the induced changes in another 
set of variables. This study represents a different approach, in 
that it utilizes computer modeling in the development of a data base 
of the relationships under investigation. 

1.1 Statement of the Problem 

The principle hypothesis of this study is: 

- Plant canopy reflectance and atmospheric modeling can 
be used to infer intrinsic scene geometry variables 
from spectral measurements. 

In dealing with this principle hypothesis two sub-hypotheses can be 
identified: 

- Wheat canopy reflectance can be predicted by computer 
modeling both as a function of low and high plant 
densities and as a function of sun angle. 

- Low wheat density and high wheat density categories 
can be inferred from model data sets. 

The principle hypothesis stipulates that a procedure can be 
developed which utilizes modeling in relating composite scene spec- 
tral measurements to the geometric make-up of the scent. The vari- 
ables used to identify intrinsic scene geometry are: 1) leaf area 

index (LAI), which is a measure of plant density; 2) leaf angle 
distribution (LAD), which is a measure of plant orientation; and 
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3) spatial dispersion, which is a measure of foliage clumping. This 
study specifically involves relating wheat plant density (i.e., LAI) 
to spectral measurements. 

The two sub- hypotheses identify specific aspects involved in the 
development of the classification procedure. The first sub-hypothe- 
sis deals with the validity of applying a canopy reflectance model to 
wheat. The second addresses the limits which can be tested for in- 
ferring scene geometry through computer modeling. 

The development of these hypotheses is an outgrowth of several 
previous investigations undertaken at Colorado State University. 

These studies can be divided into two broad categories: 1) canopy 

reflectance modeling; and 2) biomass mapping from spectral measure- 
ments. The first category is represented by R. Oliver, J. Smith 
and the author's work in developing a stochastic canopy reflectance 
model tor natural grassland (Smith and Oliver, 1972, 1974). The 
study specifically addressed the development of a computer model 
which could account for bi-directional effects of reflectance from 
a Blue Gramma canopy type. This orientation resulted in a primary 
concern for sun angle effects and canopy reflectance. 

The second category is represented by J. Tucker's, R. Pearson's, 
L. Miller's and G. Johnson's studies into the inferrence of plant 
biomass from remote sensing data (Tucker, 1973; Pearson and Miller, 
1972, 1973; Johnson, 1975). Their studies again dealt with natural 
grassland canopy types. TucKer was primarily concerned with surface 
measurements, Pearson and Miller with aircraft data, and Johnson 
dealt with satellite data. 
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The justification for merging the thrusts of these basic studies 
is the difficulty in obtaining and controlling empirical data rela- 
ting canopy spectral variations and biomass. The modeling efforts 
have the potential of deriving data sets which effectively control 
the scene variables and provide a framework for field measurement 
activities. 

The vehicle for this study was a NASA sponsored project prin- 
cipally concerned with the application of the canopy reflectance 
model in developing a wheat signature extension algorithm for sun 
angle. This was an applied project which was an extension of 
Colorado State University's earlier efforts in sun angle reflec- 
tance variations. Oliver's earlier work had resulted in a canopy 
reflectance model which indicated the potential for simulating sun 
angle dependence. The NASA effort was designed to (1) apply the 
model to wheat and (2) develop a formal approach for sun angle cor- 
rection algorithms. 

However, it was felt that the field measurement program could 
permit an initial evaluation of the capability of a modeling ap- 
proach to infer intrinsic scene parameters (e.g., LAI). If the 
results warranted further study, this initial investigation could 
serve as a guide for a more applied, in-depth experiment. 

1.2 Approach 

Figure 1 depicts a schematic of the approach used in this study. 
It consists of three primary activities: 1) field measurements; 

2) model simulations; and 3) leaf area index classification. The 
field measurements program was designed to determine the values of 



5 



FIGURE 1. THESIS APPROACH. The approach of the thesis consists of the 
general activities of Field Measurements, Model Simulation and Clas- 
sification. 
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the input parameters for the canopy reflectance model and to measure 
canopy reflectance for subsequent model calibration and evaluation. 

The model simulation activities resulted in a data base incorpor- 
ating spectral measurements as a function of leaf area index. The 
basic reflectance data was derived by executing the canopy reflec- 
tance model with a nominal input data set, with systematically 
varying LAI. In addition, an evaluation of the model's performance 
was made. An atmospheric model was employed to translate both the 
field and model derived canopy reflectances into simulated satellite 
radiance. Sun angle corrections were made to the empirical data in 
order to make it compatible with the sun angle used in the canopy 
reflectance model. 

The final activity consisted of the actual classification of 
leaf area index from spectral measurements. This process involved 
the mathematical comparison of the spectral signature of an area 
in which LAI was to be determined, to each of the modeled signa- 
tures. The leaf area index of the area was inferred to be the same 
as that of the "closest" model derived signature. Several data 
transformations were employed in order to investigate alternative 
classification procedures. 

1.3 Thesis Organization 

The origanization of this paper follows the principle activities 
outlined above. The following section identifies other research 
which is related to this study. Sections 3.0 through 5.0 are con- 
cerned with the field measurements program. These sections contain 
considerable discussion of the procedures used in assessing the leaf 



7 


angle distribution of a plant canopy. They are particularly signif- 
icant as they describe and evaluate two new techniques developed 
during this study. 

Section 6.0 is concerned with the model simulation effort. A 
brief discussion is devoted to the concepts of the two models, and 
the results are presented and interpreted. Section 7.0 discusses 
the classification algorithm, and Section 8.0 reports the results 
of the classification for leaf area index. The final section sum- 
marizes the conclusions which can be drawn from the study, and 
identifies areas of future research. 



2.0 BACKGROUND 


The activities involved in this study can be grouped into two 
broad categories: ground truth collection and extensive area map- 

ping. The first category deals with the necessity to calibrate the 
canopy reflectance model with accurate estimates of leaf area index 
(LAI) and leaf angle distribution (LAD). Existing field techniques 
were inappropriate in this application as they were either too 
tedious or yielded only approximate estimates. The second category 
is the association between target status and remote sensing data. 

The bulk of previous research in this area has been empirically 
based, and has resulted in mixed success. 

2.1 Field Procedures for Estimating LAI and LAD 

The more traditional approach to assessing LAI and LAD of an in- 
dividual plot is by direct measurement (Suits, 1972; Vanderbilt, 
1975; Oliver and Smith, 1973). The procedure involves first clip- 
ping all of the plants within a small area, then recording their 
geometric parameters. The principal measurements include plant 
height, number of leaves, length and width of each leaf, and the 
series of angles from the horizontal used to approximate the leaf's 
curvature. From these measurements the plant surface area in the 
plot can be readily determined, which, in turn, can be directly 
translated into leaf area index. The leaf angle distribution for 
the plot is estimated by calculating the normalized frequency of 
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pccurrance of each leaf inclination angle, weighted by the length of 
the leaf segment forming the angle. The geometric parameters of an 
entire field are developed by averaging a representative number of 
plot statistics. 

The direct measurement method is potentially the most accurate. 
However, it has two serious drawbacks. First, it is extremely te- 
dious and time consuming, which tends to minimize the number of 
plants used in developing plot statistics. The second short-coming, 
which is not bounded by experimentor fortitude, is ingrained in the 
non- in situ procedure of measurement. The physical removal of the 
plants can greatly affect the estimate of LAD. Wilting and loss of 
true stem orientation are obvious problems. LAI assessment is par- 
ticularly confounded by the necessary determination of discrete 
plot boundries. The principal advantage of the technique is its 
conceptual and mechanical simplicity. 

Variations of the direct measurement method include ocular es- 
timate extension and the application of a surface area meter. The 
ocular or visual method of estimation requires the visual observa- 
tion of several training plots by a technician (Pechanec and Pick- 
ford, 1937). The LAI of these sacrifice plots are then measured in 
order to refine the observer's estimations. After an adequate 
training period, the trained observer's estimations can be used to 
augment the direct measurements. This extension, however, suffers 
from human variations among estimators, is limited in spatial and 
temporal extent of the application, and does not directly character- 
ize the vegetation. 
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Another extensively used technique for determining canopy geome- 
try is the point quadrat method (Goodall, 1952; Wilson, 1973; Knight, 
1970). This method involves the calculation of LAI and LAD by de- 
termining the average number of pin contacts with the plant canopy 
for a slender, sharp-pointed rod oriented at a fixed angle with 
reference to the ground surface. The average number of contacts 
for several pin-angles are then substituted into derived linear 
equations relating LAI and the leaf inclination angle to the fre- 
quency of contact. The primary advantage of the point quadrat 
technique is its in situ nature, which minimizes plot disturbance. 

The major disadvantages are the substantial amount of time involved 
in characterizing each sample plot and sensitivity of the derived 
equations to genetic, environmental and phenol ogi cal changes. A 
more detailed discussion of the technique is presented in Section 
5.0. 

Two other field methods of LAI determination warrant brief 
consideration. They are direct measurement by a capacitance meter 
(Fletcher and Robinson, 1956) and dry weight of plant material es- 
timation (Harlan, 1976). Both procedures are driven by cfianges in 
biomass, which can be converted to LAI estimates through a double- 
sampling method described by Wilson (1963). 

t Capacitance meters measure the mass of the vegetation between 
two or more metallic probes of a specially designed capacitor which 
are inserted into the plant canopy. The meter utilizes the signif- 
icant differences between the dielectric constant of air (c'\-l), and 
that of water (e''..80). As the density of plant material increases, 
leaf water present in the vegetation increases, which yields higher 
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capacitance readings. These readings must be calibrated for varying 
species composition and environmental conditions. The primary ad- 
vantages of this method include simplicity of use and in situ measure- 
inent. The primary disadvantage is the induced experimental error 
caused ty the variations in soil water in the near-surface soil 
layers. 

The dry weight procedure attempts to develop a regression rela- 
tionship between plant surface area and dry weight. Like the capac- 
itance method, this technique is sensitive to changes in composition 
and conditions, and therefore requires frequent calibration. 

2.2 Spectral Methods for Estimating Biomass 

Spectro-optical methods for assessing vegetation biomass have 
oeen studied by several investigators. The research activities can 
be divided into two related camps. Those which concentrate on ground 
level canopy spectral reflectance (Pearson and Miller, 1973; Tucker, 
1973), and those based on elevated platform scene radiance (Maxwell, 
1975; Johnson, 1975; Miller and Pearson, 1971). These empirical 
studies involve the development of predicting equations through the 
joint consideration of spectral response and biomass of individual 
study plots. 

A fundamental assumption is that the electromagnetic energy stri- 
king a plant canopy is, in part, spectrally modified by the biomass 
present in the scene. It is further assumed that this uniqu oding 
can be decoupled from random signrl moise, and is sufficient infor- 
mation for biomass classification. This contention is closely re- 
lated to the hypothesis under study in this report. The salient 
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difference between the previous investigations and this one is in 
their methodologies. The earlier studies utilized empirical methods 
to develop predictive relationships, whereas this study employs 
model derived data to capture the information. Another important 
difference is the type of target under study. Whereas earlier 
studies dealt with natural grasslands, this study concentrates on 
a monoculture. 

The results of the ground based investigations have been success- 
ful in relating spectral reflectance of natural grasslands to esti~ 
mates of biomass. Classification of biomass from elevated platforms 
(aircraft and satellite) has been relatively less successful (Maxwell 
1975; Johnson 1975). A major reason for this is the inability of the 
researchers to effectively control the composition of their study 
plots. As the altitude of the sensor increases, tne area within 
the field of view also increases, which tends to increase the heter- 
ogeneity of the natural scene. The result is that the differences 
in responses that are presumed to be driven by variations in biomass 
alone, may be, in part, a function of species composition and canopy 
status. 

Another contributing factor to the increased difficulty of using 
elevated platforms is the enlarged bandwidth of the sensors. As 
the altitude is increased, atmospheric effects tend to rapidly at- 
tenuate and degrade the signal. In order to maintain a strong signal 
the spectral bandwidth is increased, thereby increasing the total 
energy reaching the sensor. This design has the affect of averaging 
the spectral signature, and conceals some of the potential narrow- 
band information. 



3.0 DATA COLLECTION AND REDUCTION 


5.1 General Description 

The field measurements activity presented in this report was sup- 
ported by the National Aeronautical and Space Agency (NASA), in con- 
nection with its Large Area Crop Inventory Experiment (LACIE). The 
entire field measurement program is ongoing and involves several 
research institutions. Among the principle groups participating in 
the data acquisition reported in this study were Colorado State 
University, Texas A & M University, Purdue University, the Environ- 
mental Research Institute of Michigan, Earth Observations Division 
of NASA, and the U.S.D.A. Crop Re’^orting Service. The program's 
activities were divided between an extensive empirical investiga- 
tion of the field spectral signatures for wheat and several related 
crops, and the measurement of intrinsic wheat scene reflectance 
variables . 

Two study sites were selected; one in Finney County, Kansas, 
which typifies a winter wheat region, and a second in Williams 
county, North Dakota, which is representative of spring wheat. The 
necessity for two sites is based on the dramatic differences in ap- 
pearance between the two groups of wheat, as governed by fundamental 
differences in crop development patterns and structure. In addition 
to this genetic stratification, periodic measurements were made to 
correspond to the major phonological stages in the crop's development. 
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Only a portion of this vast data set is utilized in this report. 
The data is from the Finney County site during the March 20 and 
April 23 recording periods, which correspond to the tillering and 
heading stages of winter wheat. Appendix A presents these data. 
Figures 2 and 3 are descriptive photographs of these periods. The 
striking differences between the total amount of plant cover in each 
should be readily apparent. Less obvious is the change in the leaf 
angle distribution toward a more erect canopy in April. Other sub- 
tle changes in the intrinsic parameters include a drying out of the 
soil (lighter color) and a slight decline in over-all leaf vigor 
(incipient chorolysis) in the transition from the tillering to 
heading stages. 

3.2 Field Measurement Procedures 

The field measurement procedures used in this report can be sub- 
divided into radiometric and geometric methods depending on whether 
they are involved with the estimate of optical or geometric intrin- 
sic variables. The former group •'ncludes measurements of canopy and 
soil reflectance, global and sky irradiance, and individual leaf 
transmission. The geometric procedures include an estimate of leaf 
area index (LAI) and leaf angle distribution (LAD). The format for 
discussing the data collection process adheres to the actual sequen- 
cing of activities used in the field. 

, All of the modeling directed data was collected from a single 
field. As the method for assessing plant surface area was destruc- 
tive, a new serip*; of plots had to be established for each reporting 
period. Plot selection involved the establishment of three 2' by 2* 



FIGURE 2. MARCH FIELD MEASUREMENT SITE. Rela- 
tively low plant density typlified the March 
field measurements. A pronounced "rowing" ef- 
fect between bare soil and vegetation was pres- 
ent. 



FIGURE 3. APRIL FIELD MEASUREMENT SITE. A sub- 
stantial increase in plant density, coiTipared 
to the March period, was present during the 
April field measurements. 
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plots throughout the field, which were chosen to typify the expected 
variance in the field. A rigorous statistically based sampling de- 
sign was not used because of the minimal number of plots involved. 

It was felt that an experienced agronomist could more effectively 
mentally assess the field and position the plots. A 2' by 2' con- 
figuration was chosen for two reasons: 1) it allowed for the con- 
sideration of two adjacent rows {drill spacing was approximately 
10"); and 2) it entirely bounded the field of view of the radiometer 
used to record canopy reflectance. The objective of this design and 
positioning was to identify plots which could be intensively des- 
cribed. 

3.2.1 Radiometric Measurements 

Following the establishment of the sample plots, canopy reflec- 
tance measurements were made for the remainder of the day, or until 
cloud cover became excessive. The determination as to whether atmos- 
pheric conditions were acceptable was primarily subjective, however, 
if repeated measurements of global irradiance fluctuated by more 
than 15%, it was a clear indication that radiometric measurements 
should cease. 

All of the radiometric data were collected using a LANDSAT 
radiometer commercially available through EXOTECH, INC. (Figure 4). 
This instrument is a portable, battery powered radiometer which 
uses interference filters to measure the signal in the discrete 
wavelength bands of .5-. 6, .6-. 7, .7-. 8, and .8-1. 1pm wavelength. 

The first two bands are sensitive to the green and red portion of 
visible light respectively, whereas the latter two bands record 
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responses at the slightly longer v/avelengths in the reflective infra- 
red region of the electromagnetic spectrum. These bands correspond 
to those aboard the LANDSAT-B satellite currently in orbit. The 
only modification installed on the factory unit was to convert the 
normal analog readout to a digital display in order to simplify man- 
ual recording. 

The operating procedure for determining canopy reflectance in- 
volved taking a set of readings for each wavelength band from a barium 
sulfate coated standard, then from the wheat canopy itself, "^he high- 
ly reflective and diffusing standard affords a measure of the total 
irradiance impinging on the scene, while the canopy measure identifies 
the portion of that energy which is reflected from a target. The in- 
strument's design employs a silicon cell detector positioned behind 
an interference filter and focusing optics. The ratio of the energy 
coming off the canopy to the total irradiance {outgoing/incoming) is 
defined as the reflectance of the scene for the given conditions. 

During a single "set-up" over one of the plots, four canopy re- 
flectance "observations" would be collected. These included two 
instrument positionings centered over a row, and two centered between 
rows. This multiple sampling was required as there is a separate 
aperture for each detector, which results in each having a slightly 
different portion of the canopy in its field of view. By averaging 
the four readings most of the variance due to crop rowing could be 
eliminated. A complete set of canopy reflectance measurements for a 
recording period was accomplished by rotating between the three plots, 
with a typical day consisting of twelve "set-ups" (Figure 5). 
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Also included at each set-up was a measure of the relative propor- 
tions of direct solar and diffuse sky irradiance. The procedure en- 
tailed obtaining readings of a fully exposed reflectance standard 
and then shading the standard. The unshaded measurement represents 
total global irradiance, while the shaded reading estimates the dif- 
fuse sky component. The difference between the two identifies the 
direct solar contribution. Periodically during the radiometric 
measurement period, the reflectance for bare soil in the immediate 
vicinity would be obtained in a manner similar to canopy reflectance 
determination. 

The final radiometric variable to be measured in the field was 
the actual spectral transmission of several individual leaves. An 
attempt was made to sample healthy green leaves, chlorotic (yellow- 
ing), and dead material whenever convenient. A special attachment 
to the LANDSAT radiometer was designed and constructed at Purdue 
University which enabled the direct measurement of leaf transmission 
(Figure 6). The device consists of a short cylinder (barrel), a 
sphere, and two interlocking flat disks with small slots cut in them. 
The sphere is internally coated with barium sulfate reflectance 
paint and contains a blocking baffel to prevent the passage of direct 
solar radiation, thereby, insuring only diffuse radiation at the de- 
tector. The entire apparatus fits over a single sensor port and 
must be shifted for measurements in each of the four bands. 

The sequence of operational steps begins with the alignment of 
the unit toward the sun through the use of a pinhole type sight on 
the side of the barrel assembly. Once in alignment, an unobstructed 
aperture reading is taken. The leaf is then inserted between the 



19 


two disks, so as to completely obscure the slots, and securely locked 
in place. A second reading is made which represents the proportion 
of light transmitted through the leaf. The ratio of the second read- 
ing to the first determines the faction of the total incident energy 
in a band which is transmitted through the leaf. 

Constituent transmission measurements could not be made for the 
March and April periods under study in this report. An estimate of 
this measure was made by using the green leaf transmission obtained 
during the May field session. The large proportion of green leaves 
noted during the LAI measur'^ments of both periods tends to support 
this data substitution. Also inherent in this procedure is the as- 
sumption that healthy green leaves at one crop development stage 
exhibit similar optical properties as those in another stage. The 
close agreement between model prediction and field measurement in- 
dicates that this was a reasonable assumption for the plots used in 
this study. 

Individual green leaf reflectance was approximated by setting 
these values equal to the measured transmission values. This rela- 
tionship of constituent transmission and reflectance for healthy 
green leaves being nearly equal, holds for most plant types. A 
study by Gausman (1971) shows a similar relationship in the wheat 
plants they studied. A problem in the application of this assump- 
tion was encountered for constituent measurements in band 7 (.9 - 
l.lym). The average transmission value was calculated as 54. 6T. In 
order to conserve mathematical integrity, both transmission and re- 
flectance values were set at 49.5;.. A special attachment for the 
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LANDSAT radiometer which directly measures constituent reflectance 
is currently being designed, and should be in use during the 1976 
field season. 

3,2.2 Geometric Measurements 

Field measurement of plant surface area in each plot involved a 
tedious procedure in which the one-sided surface area of all the 
living plant material in the 2' by 2' plot was determined. The ratio 
of the total one-sided surface area of all photosyntheti cal ly active 
leaves to a unit area of ground defines leaf area index. This measure 
can be easily conceptualized by imagining a tall rectangular box being 
dropped edge-wise over the canopy, without any disturbance. The total 
one-sided surface area of all living plant material (stalks or leaves) 
encased in the imaginary box, divided by the surface area of the bot- 
tom edge of the box defines leaf area index as used in this report. 

The first step in determining LAI was to remove all standing mat- 
erial within each plot. An assumption was made that equal amounts of 
plant material extended into, as out of, the plot. This simplified 
the extraction task by limiting consideration to only those plants 
whose bases were within the plot boundary. The bagged material was 
transported to an appropriate facility for surface area measurement. 

A portable Lamda Surface Area Meter (Lamda Instrument Corporation) 
was utilized for direct measurement of leaf surface area. The instru- 
ment consists of a rechargeable power supply, digital readout, and de- 
tector head (Figure 7). The detector head contains a six inch long 
bank of photo cells and a matching bank of light sources. As a leaf 
is passed between the sources and detectors, the leaf will interdict 
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FIGURE 4. LANDSAT RADIOMETER. 
The field radiometer used in 
this study is sensitive to 
four broad wavelength bands 
in the visible and near in- 
frared portions of the elec- 
tromegnetic spectrum. 



FIGURE 5. FIELD MEASUREMENTS 
OF CANOPY REFLECTANCE. Al- 
ternating radiance readings 
from the wheat canopy and a 
barium sulfate coated ref- 
erence determined canopy 
reflectance. 



FIGURE 6. LEAF TRANSMISSION 
ATTACHMENT. A special at- 
tachment to the LANDSAT 
field radiometer was con- 
structed to determine indi- 
vidual leaf transmission. 



FIGURE 7. SURFACE AREA METER. 
The surface area of plant 
mater. was determined by 
passing the material in an 
acetate carrier through the 
photo sensitive detector of 
the instrument. 
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some of the light paths, thereby indicating the total width of all 
plant material in the instrumeni; at a point in time. The total 
surface area is the summation of all the areas of the rectangular 
segments formed by the instrument's descrete sampling, as the mat- 
erial is pulled through the detector head. 

The procedure first involves separating all of the plant mater- 
ial in each field plot into live and dead components. The surface 
area of all the living material was then measured using the surface 
area meter. The ratio of the total leaf surface area to the surface 
area of the field plot identified the leaf area index of the plot. 

In most cases, the determination of L'\I consumed an entire day 
for each plot. The field measurement techniques for estimating leaf 
angle distribution proved much more expeditious. Two techniques were 
employed. The Fredholm method is based on a series of off-angle 
photographs of the wheat canopy (Figure 8). This method could only 
be used during the March 20th period, as it requires a relatively 
low plant density. The actual field procedure consists of taking 
color slides at 0, 10, 20, 30, 40 and 50 degrees off-normal; both 
with and perpendicular to the direction of crop rows. 

The Fourier technique requires several horizontal field photo- 
graphs of a thin portion of the plant canopy, silhouetted against a 
white backdrop (Figure 9). Ten inch squares were drawn on the back- 
drop to facilitate procedures involved in data reduction. A com- 
plete set of photographs for a plot consisted of three photos taken 
perpendicular to the field's rows and three aligned with the direc- 
tion of rowing. Fourier diffraction analysis, as explained later, 
was then applied to the images. 
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FIGURE 8. FREDHOLM FIELD TECHNIQUE. The field 
procedure for the Fredholm method consists of 
a series of multiple view angle ohotographs 
of the plant canopy. 



FIGURE 9. FOURIER FIELD TECHNIQUE. The Fouri 
field procedure consists of orthogonal pairs 
of silhouetted photographs. 
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3.3 Field Data Reduction 

A computer program was developed which analyzes the field radio- 
metric measurements and calculates the simple canopy reflectance, 
soil reflectance, direct to total irradiance ratio, and constituent 
transmission. This information was combined with the derived leaf 
area indices for each plot and inferred constituent reflectance to, 
yield a data base stratified by field plot and phenological stage. 

The determination of leaf angle distribution, however, required sub- 
stantial additional data reduction. The remainder of this section 
discusses the procedures involved in this process, whereas the fol- 
lowing sections describe the theoretical framework and evaluation of 
the techniques employed. 

The technique for assessing leaf angle distribution which utilizes 
multiple off-angle photographs of the plant canopy was recently de- 
veloped by R. Oliver and J. Smith (1974). This procedure solves a 
Fredholm integral to estimate LAD, and has been termed the "Fredholm 
technique." 

In general terms the Fredholm algorithm capitalizes on the readily 
observable phenomenon of increasing apparent plant density with in- 
creasing view angles. This response is primarily due to the dimin- 
ishing probability of foliage gap (seeing through the canopy to the 
soil) with views nearer to the horizon. When looking straight down 
the scene is marbled with patches of bare soil and dead material. As 

I 

the line of sight becomes more glancing, more of the live plant mater- 
ial dominates the scene. The surface area of the standing material 
projected in the direction of the camera, increases as view angle in- 
creases. The rate and pattern of this increase in projected surface 
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area serves as information for the Fredholm integral estimation of 
a canopy's leaf angle distribution. 

The data reduction procedure involved with this technique re- 
quires the determination of percent foliage cover in each off-angle 
photo. This was accomplished by overlaying a transparent dot grid 
on each print and recording the proportion of dots which did not in- 
tersect a foliage element. The probability of gap in each of the 
photographs serves as input to a computer algorithm which estimates 
the distribution of angles in the plant canopy. 

It should be apparent that the Fredholm technique is inappro- 
priate for dense canopies. At the extreme, the canopy can become 
so dense as to prohibit any foliage gaps, even when viewed from the 
vertical . 

A radically different procedure had to be developed for assessing 
LAD at crop development stages with relatively high plant densities. 
Fundamental to its methodology is the optical generation of a Fourier 
diffraction pattern, which resulted in its being termed the "Fourier 
technique". 

This method consists of two fundamental steps. The first step 
involves the determination of the distribution of angles formed by 
the foliage elements in a pair of orthogonal projections of a plant 
canopy. The two distributions are established by sampling the angu- 
lar bias of the diffraction patterns generated from the orthogonal 
projections. These diffraction patterns can be conceptualized as a 
statistical surnnary of all the orientations and arrangements in the 
original scene. The second major step involves the mathematical 
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convolution of the two planner distributions of angles into an es- 
timate of the canopy's 3-space leaf angle distribution. 

The operational steps involved in the Fourier technique's 
data reduction begins with the generation of diffraction patterns for 
each of the silhoutted photographs (Figure 10). This is accomplished 
by placing the 35mm negative of the scene into the diffractometer as 
shown in Figures 11 and 12. The resulting diffraction pattern is 
photographed (Figure 13). Experinentation proved that high contrast 
black and white film was best for all the photography involved with 
this technique. 

The angular bias of the 35mm negative of the diffraction pattern 
is then sampled with a photo cell densitometer (Figure 14). A wedge 
blocking filter is attached to the detector which limits its field of 
view to a pie shaped wedge emanating from the center of the diffrac- 
tion pattern. As the wedge filter is rotated, the intensity of 
light passing through the negative is recorded for the si'ccessive 
"pie slices". The measured distribution of angles in the diffrac- 
tion pattern acts as input to a computer program (program PROP, des- 
cribed in Appendix C) which calculates the distribution of leaf 
angles in the original scene (Figure 15). 

The final step involves the execution of another computer program 
(program CONVOL, described in Appendix C) which uses the distribution 
of leaf angles in two orthogonal projections in estimating the 3-space 
angle distribution of the canopy. 




FIGURE 10. FOURIER FIELD 
PHOTOGRAPH. The "input" in 
the generation of an optical 
diffraction pattern is a 35mm 
high contrast negative of 
silhouetted plants. 


FIGURE 12. INPUT PLACEMENT IN 
DIFFRACTOMETER. The field 
negative is positioned be- 
tween the LASER and the 
first lens. 


FIGURE 11. LASER DIFFRACTO- 
METER. The optical diffrac- 
tion pattern of the "input" 
is displayed on the screen 
at the end of the optical 
bench. 


FIGURE 13. WHEAT DIFFRACTION 
PATTERN. A diffraction pat- 
tern can be conceptualized as 
c. statistical summary of all 
the orientations and spacings 
in the "input". 
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FIGURE 14. DIFFRACTION PATTERN SAMPLING. A 
photo sensitive probe with a wedge filter is 
rotated around the diffraction pattern, while 
discrete readings are made. 





FIGURE 15. WHEAT DISTRIBUTION OF ANGLES. A 
computer program analyzes the diffraction 
pattern sampling data to infer the distri- 
bution of angles in the "input". 


4.0 THEORETICAL FRAMEWORK OF LAD FIELD TECHNIQUES 

The development of field techniques for assessing the leaf angle 
distribution of a plant canopy is a major effort of this study. The 
methods used by previous investigators yield incomplete data or are 
time consuming. This section is designed to familiarize the reader 
with the two new procedures employed. Considerable emphasis is given 
to the Fourier technique, as it was developed during this study. Min- 
imal attention is directed toward the theoretical foundation of the 
Fredholm technique, as it is aptly discussed in the report by Oliver 
and Smith (1974). A sensitivity analysis of the technique, however, 
is presented in order to identify its more important aspects. 

4.1 The Fredholm Technique 

The Fredholm technique utilizes the rate and pattern of the in- 
crease in projected surface area noted at multiple view angles of a 
plant canopy, to infer the leaf angle distribution of the canopy. A 
relatively errect canopy will exhibit a slow progression of increasing 
plant projected surface area as the view angle is increased towards 
the horizon. An Inflection point in the trend is usually observed. 

An opposite trend and pattern is noted for a flat canopy. 

The proportion of gap, PpCe^). as a function of view angle is de- 
pendent on the mean canopy projection in the direction of view, aver- 
aged over all foliage elements. Several explicit expressions of this 
functional dependence are given in the literature (Nilson, 1971). For 

example: 
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( 1 ) 


where g(6p) Is the mean canopy projection in the direction e^; and 
LAI is the leaf area index. Given a measured Pq(o). we can then in- 
vert the expression to derive g(6p). 

This mean canopy projection in direction g(0^) can then be rela- 
ted to the leaf slope distribution f( 0 g) via a Fredholm integral equa- 
tion of the first king (Oliver and Smith, 1974); g(0^) = /y^K( 0 , 0 a) 
f( 0 )de, where the kernel K(0,0g) takes a different form depending on 
whether 0=, < ~ -0., or 0. > - - A numerical solution to this 
equation, given a measured PQ{0y.), has been implemented in FORTRAN, 
and is presented in the article by Oliver and Smith. 


4.1.1 Sensitivity Analysis of the Fredholm Technique 

The sensitivity analysis procedure developed for understanding the 
Fredholm technique is unique, as it studies the effects of changes in 
an input vector on an output vector. Classical sensitivity analysis 
involves incremental perturbations to a single variable of a complex 
model, while noting the induced changes in the model's prediction. 

This relationship is usually expressed in graphical form as a percent 
change in the output verses the percent change in the input (all 
other input variables are set at nominal values). 

In the case under study, a methodology had to be derived which al- 
lows for control over the incremental changes in the input vector, 
and enables a description of the output vector changes. Figure 16 
is a general flow chart of the procedure used. The probability of 
foliage gap in multiple view angle photos (P(6AP)) forms a nineteen 
element input vector to the Fredholm algorithm. The first element 



FIGURE 16. APPROACH TO FREDHOLM SENSITIVITY ANALYSIS. Regression equations for 
both the input and output vectors constituted the basis of the approach. 
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identifies the probability of gap frwi a vertical view, with the re- 
maining elements defining off-angle views to the horizontal in incre- 
ments of five degrees. Actual physical changes in canopy leaf angle 
distributions (LAD) are characterized in unique changes in the P(GAP) 
vector. The sensitivity analysis of the Fredholm technique investi- 
gates the responsiveness of the algorithm to controlled changes in 
this input vector. 

The procedure begins with the fitting of a nominal input vector 
with a prediction equation (Figure 17). A perturbated data set is 
constructed by inducing scaler, rotational and translational shifts 
through uniform changes in each parameter of the prediction equation 
(Figure 18). The model is then run with each of the artificial in- 
put vectors, and the output noted. 

The analysis of the output involves a procedure similar to that 
of generating the perturbed data sets. The model is executed with a 
nominal input vector, and a prediction equation is fit to its output. 
A gross analysis of the model's sensitivity is obtained by plotting 
the average percent deviation from the nominal output associated with 
the percent change in the input vector. Figure 19 shows this result 
for the three types of modifications. Note that the Fredholm tech- 
nique is most sensitive to translational shifts in the P(GAP) vector, 
and nearly insensitive to scaler shifts. 

Before this analysis can be meaningful, physical factors must be 
identified with each form of input modification. An upward scaler 
shift in the P(GAP) vector implies that the probability of seeing 
soil in all of the multiple angle views is increased. This effect 
could most likely arise by an overall decrease in plant density. 
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FIGURE 17. REGRESSION FIT OF INPUT VECTOR. 
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FIGURE 18. INPUT VECTOR MODIFICATIONS 
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as measured by the leaf area index. Clockwise rotation of the input 
vector results in an increase in the probability of soil at the more 
vertical view angles, with a simultaneous decrease at the more hor- 
izontal views. This effect is most likely the result of the canopy's 
LAD becoming more vertically biased. Translational shifts represent 
relative plant dispersion and spacing which are driven by complex 
interactions of both LAI and LAD. It can be generalized, therefore, 
that the Fredholm technique is most sensitive to changes in canopy 
dispersion and LAD, and relatively insensitive to LAI, until plant 
cover becomes extremely dense. 

Figures 20 through 22 show the type and magnitude of the induced 
effect on the model's output. The "a" and "b" curves identify scaler 
and rotational shifts in the output vector (LAD), respectively. In 
general, it can be stated that all three types of input vector mod- 
ifications (scale, rotation, and translation), manifest themselves 
primarily in "a” factor shifts of the predicted LAD vector. As this 
vector is renormalized in the actual Fredholm algorithm, a purely 
scaler shift in LAD is not possible. However, an increase in the "a" 
factor does indicate a bias towards smaller inclination angles and a 
flattening of the slope of the LAD vector. 

The "b" parameter best describes the relative form of the leaf 
angle distribution. A large positive "b" indicates a steep slope, 
which in turn, implies an erect canopy. In reviewing Figures 20 
through 22, it is apparent that translational and rotational shifts 
dominate the determination of the slope of the output vector. Spe- 
cifically, counter-clockwise rotation of the P(GAP) vector results 
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FIGURE 19. AVERAGE SENSITIVITY. Translational 
shifts in the P{6AP) input vector had the 
greatest effect on the LAD output vector. 



FIGURE 20. SENSITIVITY TO SCALER SHIFTS. 
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FIGURE 21. SENSITIVITY TO ROTATIONAL SHIFTS. 




FIGURE 22. SENSITIVITY TO TRANSLATIONAL SHIFTS. 
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in an increase in the "b" parameter. The movement of the inflection 
point to the lower view angles generates a similar response. 

In summary, the Fredholm model's prediction of leaf angle distri- 
bution is most responsive to rotational and translational shifts in 
the P(GAP) input vector. Clockwise rotation of the P(GAP) vector 
is principally the result of a flattening of the foliage's inclina- 
tion angles, and results in the LAD' s bias towards the smaller 
angles. Movement of the P(GAP) vector's characteristic inflection 
point to the left implies a complex condition (dependent on both 
LAI and LAD), where the projected surfaces of the foliage elements 
tend to become rapidly obscured at smaller view angles. This condi- 
tion denotes LAD'S bias towards the larger inclination angles. 

4.2 The Fourier Technique 

The Fourier technique involves the mathematical convolution of the 
distributions of angles for two or more orthogonal horizontal projec- 
tions of a plant canopy to infer the three-space leaf angle distri- 
bution. In discussing this method, the physical process generating 
diffraction patterns is first presented, followed by discussions on 
the mechanics of sampling the patterns and the theoretical basis of 
convoluting the orthogonal distributions of angles. 

4.2.1 Fundamentals of Optical Diffraction Patterns 

Diffraction patterns result from the diffraction of light rays 
and their subsequent unique interference patterns. Diffraction re- 
fers to the bending of waves around an obstacle. Interference is an 
effect that occurs when two 'aves of equal frequency are superimposed. 
If, at the point of meeting, two waves are in phase (vibrating in 
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unison, and the crest of one coinciding with the crest of another), 
they combine to form a new wave of the same frequency and a larger 
amplitude equal to the sum of the amplitudes of the original waves 
(Figure 23). This effect is termed constructive interference. 

If two waves meet out of phase (crest of one coinciding with a 
trough of the other), the result is a wave whose amplitude is the 
difference of the original amplitudes. The counterpart process is 
termed destructive interference. If the original waves have equal 
amplitudes and meet out of phase, they can be conceptualized as 
completely destroying each other; leaving no wave at all. Purely 
destructive interference, therefore, produces a dark spot, while 
constructive interference produces a bright spot. Partial construc- 
tive or destructive interference results whenever the waves have an 
intermediate phase relationship. 

The process generating diffraction patterns involves optical 
diffraction and interference (Dobrin, 1968; Goodman, 1968), The 
diffraction of light as it passes edges, uniquely changes the phasing 
of light waves, thereby altering the interference relationships form- 
ing the diffraction pattern. As the relationships uf the edges 
change, a corresponding change is induced in the interference pattern. 
Figure 24 schematically shows this fundamental process. 

Stable diffraction patterns can only be generated if the phasing 
of the light rays forming them are not constantly changing. This 
stable condition is termed coherent. Most light is non-coherent, as 
the radiation from different atoms, acting as light sources, are con- 
stantly changing their relative phase relationships. An ordinary 
light source radiates light of mixed wavelengths in an out-of-step 
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FIGURE 23, CONSTRUCTIVE/DESTRUCTIVE INTERFERENCE. 
Interference describes the Interaction between 
two waves when they met. 
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FIGURE 24, GENERATION OF DIFFRACTION PAHERNS. 
The display of the unique Interference pattern 
on a screen Is termed a diffraction pattern. 
Note the difference In diffraction pattern 
spacing and slit spacing between the two 
schematics. 




40 


(non-coherent) and random manner. In contrast* LASERS emit mono- 
chromatic, coherent light waves. The single slit in Figure 24 is 
required as it converts a normal light source into a nearly coherent 
source. In current applications, a LASER is used as the light source 
in diffraction analysis, because of its coherent properties and 
strong illuminance. 

Figure 25 shows the process of generating diffraction patterns 
with the aid of a lens. This configuration allows for the concen- 
tration of the energy in the diffraction pattern, and results in 
much stronger and more discernable patterns. 

Figure 26 outlines the fundamental relationships between the in- 
put image and the diffraction patterns which can be analyzed for 
spatial information (Pincus, 1969; Pincus and Dobrin, 1966), The 
types of information that can be uncovered include: 

1, Orientation of elements 

2, Spacing of elements 

3, Location of elements 

Element direction information is contained in the orientation of the 
diffraction "dots", which are perpendicular to the linear elements 
producing the dots. The spacing Intelligence is inferred by the 
spacing of the dots, with the distance between the dots varying in- 
versely with the true spatial frequency in the input. The tech- 
niques used in determining element location involve the related 
field of Holography and operate on the retention of the true phasing 
of two wave fronts. The work reported in this paper capitalizes on 
the ability to infer element orientation from diffraction patterns. 
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FIGURE 25. LENS CONTROLLED DIFFRACTION PATTERNS. 
The introduction of a lens between the "input" 
and screen increases the intensity of the dif- 
fraction pattern. 



FIGURE 26. DIRECTIONAL INFORMATION IN DIFFRACTION 
PATTERNS. The orientation of the diffraction 
pattern is perpendicular to the linear segments 
itt the "input". 
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In addition to the principle effects driving the orientation, 
spacing and location information, secondary effects are induced 
through multiple order considerations, interactions between ele* 
ments, and irregular element spacing and orientation. The simul- 
taneous interactioi of both the primary and secondary factors 
results in clouds of points rather than descrete dots. These 
clouds can be thought of as a statistical summary of all the spa- 
tial information contained in a complex scene. 

4.2.2 Sampling and Interpreting Orientation 

Figure 27 shows the components of a si.nple optical diffracto- 
meter used to generate diffraction patterns. The LASER acts as a 
monochromatic, coherent light source which illuminates a reduced 
transparency (input) of a field scene (in this study the scene con- j 

sis ted of a horizontal picture of silhouetted wheat). The input 
functions as a diffraction grating with unknown spatial properties 
and bends the waves in a unique way. The altered light rays are 
passed through an objective lens which focuses the diffraction pat- 
tern onto a plano-convex lens. This final lens enlarges the pattern 
and displays it onto a rear projection screen. The diffraction pat- 
tern can then be photographed and retained for later analysis. 

Sampling the diffraction pattern for spacing and orientation in- 
formation about the original scene involves recording the portion of 
the pattern in appropriate segments. Spacing intelligence, as ear- 
lier noted, is contained in the distance between diffraction dots 
and the origin of the pattern. By placing consecutively larger cir- 
cular bands around the origin and measuring the relative amounts of f 





the pattern contained in each, an estimate of the distribution of 
spacings can be obtained (Figure 28). Similar to ring sampling for 
spacing bias, is v/edge sampling, which involves rotating a "pie slice" 
shaped filter about the origin. A distribution of angles in the input 
is inferred by the relative proportions of the diffraction pattern 
contained in each slice. Simultaneous sampling with a ring and wedge 
filters results in an estimation of the proportion of lines with a 
particular spacing that are oriented in a specific direction. 

Figure 29 shows the wedge sampling procedure used in this study. 
The high contrast negative of the diffraction pattern is placed on a 
light table and covered with a small circular frame. A filter de- 
signed to fit into the frame, was construe vad from a disk with a 
three degree slice removed. The relative intensity of light passing 
through the diffraction pattern segment defined by the wedge filter, 
is measured by a digital photometer. As the filter and attached de- 
tector are advanced equal distances (corresponding to equal degree 
rotations), discrete readings of the proportion of angles are re- 
corded. 

The tonal mapping of the diffraction pattern by the photographic 
negative is reversed. Thi ' efore, a high reading by the photometer 
indicates a minimal portii/n of the pattern, while a low reading iden- 
tifies the presence of a large portion. Program PROP (Appendix C) 
analyzes this information and develops the distribution of angles 
contained in the input scene. 

In the application of spatial analysis to determining the leaf 
angle distribution of an entire canopy, it is necessary to convolute 
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FIGURE 27. DIFFRACTOMETER COMPONENTS. The com- 
ponents of a laser diffractometer consists of 
a coherent, monochromatic light source, "input", 
focusing optics, and display screen. 
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FIGURE 28. SAMPLING FILTERS. Ring and wedge 
filters are used to sample spacing and orien- 
tation of the diffraction patterns, respec- 
tively. 
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the two o»'thogona'l probability distributions inferred from the dif- 
fraction patterns. Program CONVOL (Appendix C) is used to execute 
this process. The theoretical foundation of CONVOL's algorithm is 
presented below. 


I 


Given; distribution of angles in XZ view 
distribution of angles in YZ view 
Find : distribution of angles in XYZ space 

Step 1 develop a matrix identifying the three-space angles 

associated with all pairwise combinations of descrete 
angles (e.g., 5 degree intervals between 0-90 degrees) 
which are possible in the orthogonal views. 

Subproblem : (refer to Figure 30) 

Derive a mapping function from e in orthogonal views 
to 0 in three- space. 

Given: a line in three-space with and 

assume lowest point is at origin (X=Y=Z=0) 

Solution; Zg = a + bX 2 
— 0 tan 
and if we chose X 2 = 1. 

Z 2 = 0 + tan 
= tan 0 JJ 2 

in plane YZ, 

Zg = a + bY 2 
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FIGURE 29. SAMPLING DEVICE. The components of 
the sampling device used in this study consis 
of a 3 degree wedge pass filter, a photo sen- 
sitive detector, and a digital readout. 
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FIGURE 30. LINE IN 3-SPACE 
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substituting* 

tan 0JJ2 = (0 tan 

'f 2 ~ tan Sj^^/tan ®y2 

therefore, two points on the line are known in 
three-space; 


^2 ■ ^2 " ®YZ* ^2 “ ®XZ 


using the Pathagorean therom. 


c 


^/7T7 




= V (1-0)^ t 


tan e 
tan 0 


u, 

YZ 


•• 0 


V 


1 + 


tan 9^ 2 


tan 0 


YZ 


using the tangent relationship (opposite over adjacent) 
and substituting, 

tan 0^ = (Z 2 “Zj)/c 

= (tan 0 JJ 2 " 0)/c 
= tan 


or, 


0j = tan 


(tan 0 ^^ / 


vr: 


tan 0 
tan 6 


n 

YZ 


) ( 2 ) 
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This equation can be used to reduce any pair of 
orthogonal angles describing a line in three- 
space to the three-space inclination angle of 
the line, as required in step 1. It is used by 
program THETA {Appendix C) in developing the 
data matrix used by program CONVOL. 


Step 2 The probability of any ®y2 joint 

probability of 3*1^1 ®yZ’ 

P(6xz» * ^^®YZ) 

The diffraction analysis of the orthogonal field 

photos yields the distribution of angles in the 

XZ and YZ planes. 

Step 3 A unique may occur by several planer combinations, 
therefore, 

P{6i) = E (PCe^x,) * P(eyzj)) (3) 

where s indicates summing of all of the joing prob- 
abilities of all the pairwise combinations forming 

“i- 

4.2.3 ‘ Qualitative Analysis of the Fourier Technique 

This complex procedure was qualitatively evaluated in several 
ways. A computer prograr, was designed which calculates the coor- 
dinates on orthogonal projections of lines in three-dimensional 
space. Figure 31 is a print of one of the microfilm sets of the 
two projections for a series of lines slanting at 45 degrees. It 



X VIEW 




FIGURE 31. UNIFORM ABSTRACT PLANT CANOPY EVALUATION. A computer program was developed 
which generated microfilm plots of the orthogonal perspectives of a series of lines 
with a given spatial configuration. 
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may help the interpretation of these projections by imagining then 
to be transparent adjacent sides of a cube surrounding a series of 
sticks. In the X view, the lines form about a 60 degree angle, while 
in the Y view a 120 degree angle is indicated. The convoluted dis- 
tribution of angles correctly indicates the true 46 degree three- 
space angle. Several cases were made for various angles and the 
procedure tracked them exceptionally well. 

Figure 32 shows a typical example of the same program, yet the 
positioning and angles of the “sticks" were randomly assigned in 
accordance with a pre-determined distribution of angles. It v/as an- 
ticipated that the convoluted distribution of angles would closely 
correspond to the defined distribution driving the program. Again 
promising results were noted. However, as more lines were put into 
the scene the diffraction patterns became less distinct. This is 
due to the multiple interactions between lines. In light of these 
observations it was determined that a conserted effort should be 
made to reduce the number of silhouetted plant elements in a field 
photograph. 

Figure 33 illustrates another qualitative evaluation method that 
was utilized. It consisted of taking orthogonal photographs of an 
abstract plant canopy constructed of tinker toys. All of the branches 
form 45 degree angles, with stalks perpendicular to the ground. The 
tracking of these canopy angles was excellent. 

An actual field evaluation for both the Fredholm and Fourier tech- 
niques was performed and is presented in the next section. The re- 
sults of this analysis showed the Fourier technique to be an extremely 
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FIGURE 32. RANDOM Ai ■<AU PLANT CANOPY EVALUATION. The microfilm plots of orthogonal 
views analyzed by the Fourier technique were developed by stochastic computing from 
user defined population statistics of spatial orientation. 









FIGURE 33. TINKER TOY ABSTRACT PLANT CANOPY EVALUATION. Orthogonal photographs of 
"Tinker Toy" plant canopies, with all branches set at 45 degrees inclination angle, 
were analyzed. 






accurate descriptor of leaf angle distribution. The Fredholm method 
proved to be less accurate, however, it still offers a rapid tech- 
nique for first order estimation of leaf angle distributions. 


5.0 BENCHMARK COMPARISON OF LAD FIELD TECHNIQUES 


A field evaluation was conducted which compared the Fredholm and 
Fourier methods for measuring leaf angle distributions used in this 
study. In addition, the Point Quadrat technique, extensively used in 
the agricultural sciences, and orthogonal tracings of leaves were 
employed. 

An evaluation plot was selected such that it had a low leaf area 
index (.98) and its plants were uniformly dispersed (Figure 34). The 
plant type was Western Wheatgrass (Agropyron Smithii), which is 
characterized as having a conical leaf angle distribution (Oliver and 
Smith, 1973). The plot dimensions were 16 by 18 inches. 

5.1 Description of Techniques 

5.1.1 Point Quadrat Technique 

It has been shown that the mean foliage angle can be calculated 
from the number of contacts made by point quadrats passed vertically and 
horizontally through a plant canopy (Philip, 1965; Wilson, 1959). In 
practice the error associated with this method rarely exceeds 10%. The 
technique is in situ , however, appreciable localized trampling is 
induced around the field plot. The time required for a single angle 
determination is about 18 man-hours (Knight, 1970). This method is 
most commonly used to characterize foliage geometry, however, it only 
estimates the mean inclination angle rather than a distribution of angles. 



The field procedure involves the calculation of the average num- 
ber of contacts a long slender pin makes with the vegetation during 
a pass through the plant canopy (Figure 35). The length of the pass, 
for both the horizontal and vertical transects, is dictated by the 
height of the canopy. Several hundred passes are made from both 
directions, with the averages for each being multiplied by empiri- 
cally obtained coefficients to determine the mean foliage angle. 

5.1.2 Orthogonal Tracing Technique 

The distribution of foliage angles for an individual plant can be 
accurately determined by analyzing two orthogonal photos of the plant 
(Oliver and Smith, 1974). The distribution of angles for an entire 
plot is statistically determined by averaging the distribution of 
several representative plants. This technique has many of the lim- 
iting features associated with the point quadrat method. It is slow, 
tedious and destructive. However, it is a direct method which makes 
it useful for comparing the other methods. 

With this procedure, individual plants are clipped from a field 
plot and the silhouetted profiles are photographed from two ortho- 
gonal directions (Figure 36). The photographs are then digitized 
by placin: a transparent grin over the photographs and recording the 
two-dimensional coordinates of straight line segments along the pro- 
files. The profiles are plotted on microfilm (Figure 38) using the 
digitized data in order to verify th'’ digitization. A computer 
program was developed for this study which determines the three- 
dimensional coordinates of the foliage elements from the two sets of 
orthogonal data, and calculates the average foliage inclination 



56 



FIGURE 34. EVALUATION PLOT, 
The evaluation plot was a 
relatively low density 
area of western wheat- 
grass. 



FIGURE 35. POINT QUADRAT 
FIELD TECHNIQUE. The field 
procedure consists of noting 
the number of contacts a 
long slender pin makes with 
foliage elements throughout 
numerous passes through the 
canopy at specific angles. 



FIGURE 36. ORTHOGONAL TRACING 
FIELD TECHNIQUE. Orthogonal 
photographs were made of in- 
dividual silhoutted plants. 
These same photographs were 
used for the Fourier analy- 
sis. 


FIGURE 37. FREDHOLM FIELD 
PROCEDURE. Multiple view 
angle photographs were 
made of the evaluation 
plot. 


t 
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FIGURE 38. ORTHOGONAL TRACING DATA REDUCTION. The orthogonal 
field photos of individual plants were digitized, with the 
two-space coordinates serving as input to a computer routine 
which calculated the three-space distribution of angles. 
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FIGURE 39. FREDHOLM DATA REDUCTION. The percent cover in each of 
the multiple view angle photographs was determined by dot grid 
analysis. 



angle by direct computation (program, ORTHOG, Appendix C). The dis- 
tribution for the entire plot is calculated by the weighted averaging 
of the individual plant distributions based on the size of the plant. 

5.1.3 Fredholm and Fourier Techniques 

The field procedures and theoretical framevvork for the Fredholm 
and Fourier methods are presented in the previous section. Figures 
37 and 39 show examples of the field photographs and reduced data. 

For this evaluation the tracings of the individual plant photographs 
used in the Orthogonal Tracings technique were used as input in gen- 
erating the Fourier diffraction patterns (Figure 40). 

The relative merit of the Fredholm method is its ease of data 
collection with a minimum of canopy disturbance. The Fourier tech- 
nique is much slower and more tedious than the Fredholm method, yet 
it is still relatively easy and rapid when compared to either the 
Point Quadrat or Orthogonal Tracing techniques. 

5.2 Description of Evaluation Procedure 

Two fundamentally different evaluation procedures were used. The 
first involved the comparison of the four basic techniques' predic- 
tions of the leaf angle distribution for a common field plot. The 
normal steps for the Orthogonal, Fourier, and Fredholm methods are 
schematically shown in Figure 41. The Point Quadrat technique merely 
required the simple solution of a linear equation predicated on the 
field estimation of the number of needle contacts. 

The second evaluation approach is more abstract in nature, and 
addresses the validity of the convoluting algorithm used in the 
Fourier technique. Figure 42 identifies the major steps in this 
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FIGURE 40. FOURIER DATA REDUCTION. Diffraction patterns were 
generated and analyzed of each of the orthogonal photos of the 
individual plants. 
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FIGURE 41. TECHNIQUES EVALUATION APPROACH. Each of the three 
field techniques under study were executed on a common plot, 
and compared. 
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ORTHOGONAL FOURIER 



FIGURE 42. FOURIER METHODOLOGY EVALUATION APPROACH. The convo- 
luting algorithm and orthogonal pairing assumption were evalu- 
ated. 
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process. In the left hand flows the digitized photos are used to cal- 
culate a relatively exact distribution of angles in each orthogonal 
view. The orthogonal pairs are then convoluted and finally all three- 
space distributions are averaged to yield the predicted canopy leaf 
angle distribution. The ability of this abstract approach to accur- 
ately correspond to the normal Orthogonal method is a measure of the 

validity of the convoluting algorithm. 

The right hand flows identify a process designed to evaluate the 
necessity of having truely orthogonal pairs of silhouetted field 
photos. If the canopy does not have an azimuthal bias, then there 
would be no physical reasons for differentiating between photo per- 
spectives, other than statistical sampling requirements. The ex- 
treme right hand column incorporates this assumption by determining 
the average distribution of all the planner projections, and then 
convoluting this average distribution on itself. A slightly differ- 
ent format (identified by the dotted lines) convolutes all the pair- 
wise combinations of planner distributions, and averages the result- 
ing three-space distributions. The essence of these reviews is to 
determine whether the Fourier technique is truely dependent on or- 
thogonal field photos. 

6.3 Comparison of Results 

The results of the separate procedures are summarized in Figures 
43 through 49 and Table 1. Figures 43 through 46 correspond to the 
evaluation of the normal procedures, while Figures 47 through 49 
address the convolution assumption of orthogonal views. The results 
are graphically presented to facilitate a qualitative assessment of 
each technique’s precision. 
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Table 1. COMPARISON OF FIELD TECHNIQUES FOR ASSESSING LAD. 



INCLI- 

NATION 

ANGLE 

ORTHO- 

GONAL 

TRACING 

FOURIER 

FREDHOLM 

ORTHO- 

GONAL 

CONVOLUTE 

FOURIER 

AVERAGE 

FOURIER 

PAIRWISE 

2.5° 

0. 

0. 

0. 

0. 

0. 

0. 

7.5 

0. 

0. 

0. 

0. 

0. 

0. 

12.5 

0. 

0. 

0. 

0. 

0. 

0. 

17.5 

0. 

0. 

0. 

0. 

0. 

0. 

22.5 

0. 

0. 

0.002 

0. 

0. 

0.005 

27.5 

0. 

0.001 

0.012 

0. 

0.001 

0.005 

32.5 

0.016 

0.005 

0.024 

0.018 

0.006 

0.015 

37.5 

0.026 

0.012 

0.038 

0.025 

0.015 

0.026 

42.5 

0.006 

0.027 

0.051 

0.002 

0.033 

0.055 

47.5 

0.022 

0.046 

0.063 

0.040 

0.050 

0.083 

52.5 

0.114 

0.068 

0.073 

0.058 

0.069 

0.099 

57.5 

0.056 

0.110 

0.083 

0.158 

0.105 

0.132 

62.5 

0.112 

0.121 

0.091 

0.141 

0.111 

0.119 

67.5 

0.139 

0.149 

0.100 

0.132 

0.139 

0.142 

72.5 

0.196 

0.131 

0.106 

0.129 

0.126 

0.117 

77.5 

0.158 

0.111 

0.112 

0.159 

0.112 

0.114 

82.5 

0.088 

0.130 

0.120 

0.080 

0.139 

0.069 

87.5 

0.066 

0.088 

0.126 

0.057 

0.093 

0.023 

MEAN 

FOLIAGE o 

ANGLE 68.0 

68.0° 

66.0° 

66.7° 

67.9° 

62.8° 


Point Quadrat Mean Inclination Angle = 72.5 
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FIGURE 43. ORTHOGONAL TRACING TECHNIQUE RESULTS. 
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FIGURE 44. FREDHOLM TECHNIQUE RESULTS. 
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FIGURE 45. FOURIER TECHNIQUE RESULTS. 



UAf »K|Tm K%i.l 

FIGURE 46. FREDriOLM/FOURIER DEVIATIONS. 




The mean foliage inclination angle calculated by the Point Quadrat 
method was 72.5 degrees. In comparing this prediction with those cal- 
culated by the other three techniques it appears to be an over-esti- 
mate (Table 1). This condition could easily be a result of this 
author's inexperience with its field procedure. Also contributing 
to the implied error could be that the empirically derived regres- 
sion equation employed was developed for general grassland canopies 
(Knight, 1970) and not specifically for Western Wheatgrass. 

Figure 46 compares the effectiveness of the Frednolm and Fourier 
methods to track the cummulative distribution of angles predicted by 
the Orthogonal Tracing technique. It is readily apparent that the 
Fourier procedure compares more favorably. Less obvious is the com- 
mon pattern of deviations. Both techniques tend to over-estimate 
the probability of inclination angles less than 80 degrees. This 
result is most likely due to the averaging approach both tech- 
niques which fail to respond to the sharp fluctuations of the Ortho- 
gonal method. This same effect is apparent in the predictions made 
by convoluting the actual planner distribution of angles calculated 
by the Orthogonal technique (Figure 47). 

The comparison of pairwise and average approaches of convolution 
(Figures 48 and 49) shows that strict utilization of truely orthogonal 
pairs is not necessary in these circumstances. In this evaluation, 
the convolution of the average distribution of angles for all of the 
planner views did an excellent job in tracking the orthogonal method. 
However, a rowed crop, such as wheat, may require orthogonal pairing, 
and further evaluation is warranted. 
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FIGURE 47. CONVOLUTED ORTHOGONAL TRACING RESULTS. 
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FIGURE 48. CONVOLUTED AVERAGE DISTRIBUTION OF 

ANGLES 
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FIGURE 49. CONVOLUTED PAIRWISE COMBINATIONS OF 
DISTRIBUTIONS OF ANGLES RESULTS. 



6.0 MODEL SIMULATION 


Fundamental to the hypothesis under investigation is the ability 
to adequately model the complex system of interactions between solar 
radiation and a plant canopy. In addition, the interactions between 
this radiation and the atmosphere must be linked to the canopy model 
in order to fully describe the physical process. The models adapted 
to the investigation are briefly discussed in the following subsec- 
tion. The description of these models is merely intended to famil- 
iarize the reader with their basic approaches. The cited literature 
contains several detailed reports which should be consulted for a 
more thorough understanding. 

The process used in discriminating geometric variables of a plant 
canopy principally involves the statistical comparison of a sensor 
measured response to a series of model derived responses. In devel- 
oping these model responses, the geometric variable of concern is 
allowed to systematically vary, while all of the other variables re- 
main fixed at t eir best estimates. The specific procedures used in 
simulating the surface canopy reflectance and satellite radiance arc 
presented in the last two subsections. 

6.1 Canopy Reflectance and Atmospheric Models 

The canopy reflectance model used in this study is Colorado State 
University's Solar Radiation Vegetation Canopy (SRVC) Model (Oliver 
and Smith, 1973, 1974). This model differs from other plant canopy 



models (Suits, 1972; Colwell, 1974; Allen and Richardson, 1968; 

Kubelka and Munk, 1931) in that the driving variables and canopy 
structure are based on probability distributions. This orienta- 
tion results in an estimation of both the mean canopy reflectance 
and its covariance matrix. The deterministic approach of the other 
models yields estimates of mean responses, without any direct in- 
ference to the central tendency. 

The model's input parameters can be divided into two principle 
classes: 1) environmental factors; and 2) intrinsic scene character- 

istics. The environmental factors include sun position, diffuse and 
direct irradiance, and sensor view angle. Leaf area index, leaf 
angle distribution, and spatial dispersion of foliage elements des- 
cribe a plant canopy's geometric characteristics. The canopy's 
radiometric input parameters include soil reflectance and individual 

leaf reflectance and transmission. 

The methodology of the model involves the mathematical tracking 
of a photon of light as it interacts with a plant canopy. The en- 
vironmental factors determine the spectral composition and angular 
dependency of the total scene irradiance. Intrinsic scene parameters 
are used to develop the probability distributions which stochastic- 
ally determine the interaction of the solar radiation and the plant 
canopy, and to calculate the spectral modifications and redirections 
resulting from the interactions. The model's operation begins with 
an instantaneous burst of hemispherical irradiance. This pulse is 
then charted in its multiple interactions with the canopy until all 
of it has escaped back into the sky or has been absorbed by the 
canopy and background. 



72 


The atmospheric radiative transfer model used in this study is 
the Environmental Research Institute of Michigan's Model (Turner, 
1973). The approach of this model is founded on the solution of a 
series of explicit mathematical expressions. The basic relationship 
describes the total spectral radiance at a sensor as, 


M * * T) + N 

target path 

where, N = total radiance 


(4) 


^target 


= target spectral radiance at the surface 


N 


T = spectral transmittance between the sen- 
sor and the target 
path ” spectral path radiance 


All three of these quantities depend upon the condition of the atmos- 
phere. The actual state of the atmosphere at any location and at any 
time is approximated by atmospheric variables and environmental fac- 
tors. Among the atmospheric variables are the scattering and absorb- 
ing properties of gases and particulates that exist in the atmosphere. 
Environmental factors include sensor attitude and view angle, sun an- 
gle, and canopy reflectance. 


6.2 Siii.ulation for Surface Canopy Reflectance 

The complete procedure for estimating leaf area index from spec- 
tral measurements through the use of modeling contains four major ac- 
tivities: 1) the identification of model input parameters through a 
field measurements program; 2) the construction of a model generated 
data set which tracts the induced changes in canopy spectral reflec- 
tance arising from variations in LAI; 3) the translation of surface 
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reflectance data into estimated radiance values at an airborne sen- 
sor ttirough atmospheric modeling techniques; and 4) the statistical 
comparison of measured and model derived spectral signatures to yield 
an estimate of the scene's LAI. The first activity was outlined in 
Section 3.0, and the fourth is presented in Section 7.0. This sec- 
tion is concerned with the middle two activities. 

The canopy model simulation for LAI determination utilized the 
March and April field measurements, which correspond to the tiller- 
ing and jointing stages of wheat. These stages were selected -js 
they contain the extremes in LAI and are similar in canopy make-up, 
both are characterized by a green canopy prior to seed head devel- 
opment, with their primary difference being a dramatic increase in 
LAI at the jointing stage. 

Figure 50 shows a comparison of field and model estimates of 
canopy reflectance at the two phenological stages under study. In 
both cases, the model prediction was made from a nominal input data 
set representing the best estimates of field conditions. The empir- 
ical signature is an average of the direct measurements made from 
the three intensive field plots. Additional comparisons were made 
for thirty sun angles obtained from canopy reflectance measurements 
taken between 1030 and 1630 hours on March 20th, and 1000 and 1800 
hours on April 23rd. A complete set of the field and model data from 
this effort is reported in a paper by Smith, Berry and Heimes (1975). 
In general, it can be noted that the model was successful in track- 
ing the empirical data during the jointing (April) stage, while it 
was somevjhat less accurate during the tillering stage (March). 
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MARCH . 1030 hours 



MSS BAND 
4 


5 


6 


7 


MODEL 

.050 

+.007 

.036 

.005 

.262 

.004 

.359 

.058 



MSS BAND 


MSS BAND 

4 


5 


6 


7 


MODEL 

,053 

±.002 

.037 

.001 

.281 

.030 

.390 

.014 


FIELD 

.076 

+.004 

.075 

.005 

.264 

.031 

.378 

.075 


FIELD 

.040 

+.C05 

.025 

.005 

.259 

.031 

.391 

.041 


FIGURE 50; MODEL/FIELD CANOPY REFLECTANCE COMPARISONS 
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In both cases > however, the model conformed to the trajectory of 
the empirical signature. This point is pertinent to the usage of 
data transformations in the statistical comparison process and will 
be expanded in Section 5.0. 

In simulating LAI effects, all of the variables of the canopy 
model were fixed in accordance with their field estimates, while 
the leaf area index was varied from 0.5 to 5.0 for the March simu- 
lation, and 3.5 to 8.0 for the April data. A single model execu- 
tion consisted of ten samples, each comprised of five trials. A 
.5 step in LAI was used for each execution in simulating both stages. 
The result of this process was the construction of a Simulated Can- 
opy Reflectance data set containing one hundred row entries, and 
four column entries corresponding to the four LANDSAT spectral 
bands. The row entries consist of ten model predicted spectral sig- 
natures for each of ten LAI conditions. 

Figures 51 and 52 graphically portray the relationship between 
scene reflectance and LAI for both th^ field measurements and model 
predictions in March and April. Table 2 tabularly summarizes the 
data. The model data for March identifies a positive relationship 
for MSS bands 4, 6 and 7 which appears to plateau at LAI's above 
3.5. Band 5 displays a less prominent negative relationship which 
also approaches an asymptote at about 3.5 LAI. These general trends 
favorably agree with published studies of a closely related factor, 
canopy biomass (Tucker, 1973). The empirical data for these same 
periods tend to agree with the model data, with the exception of the 


' f- 


lowest LAI in March. 
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Table 2. MODEL/FIELD CANOPY REFLECTANCE AS A FUNCTION OF LAI 


LAI MSS 4 MSS 5 MSS 6 MSS 7 MSS 7/5 


MARCH (MODEL) 


APRIL (MODEL) 


MARCH (EMPIRICAL) 


APRIL (EMPIRICAL) 


.5 

.043 

.045 

1.0 

.052 

.045 

1.5 

.047 

.037 

2.0 

.053 

.042 

2.5 

.049 

.035 

3.0 

.053 

.038 

3.5 

.050 

.035 

4.0 

.051 

.035 

4.5 

.052 

.036 

5.0 

.0496 

.034 


3.5 

.050 

.035 

4.0 

.051 

.035 

4.5 

.052 

.036 

5.0 

.049 

.035 

5.5 

.051 

.036 

6.0 

.051 

.035 

6.5 

,050 

.035 

7.0 

.049 

.034 

7.5 

.049 

.034 

8.0 

.049 

.034 


1.31 

.071 

.077 

2.07 

.076 

.078 

4.06 

.080 

.071 

5.13 

.041 

.027 

5.36 

.039 

.025 

6.15 

.040 

.025 


.149 

.200 

4.45 

.239 

.326 

7.24 

.233 

.320 

8.67 

.264 

.361 

8.63 

.257 

.353 

10.02 

.279 

.380 

9.98 

.269 

.373 

9.99 

.272 

.378 

10.29 

.277 

.385 

10.24 

.266 

.371 

10.26 


.269 

.373 

10.62 

.272 

.373 

10.68 

.277 

.385 

10.62 

.266 

.371 

10.76 

.273 

.379 

10.64 

.271 

.375 

10.65 

.269 

.373 

10.65 

.264 

.365 

10.60 

.264 

.365 

10.63 

.264 

.365 

10.63 

.298 

.475 

6.17 

.235 

.320 

4.10 

.260 

.340 

4.79 

.262 

.381 

14.11 

.266 

.401 

16.04 

.262 

.391 

15.64 











REFLECTANCE REFLECTANCE 



FIGURE 52. REFLECTANCE AS A FUNCTION OF LAI (APRIL) 
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Figure 53 shov/s that the reflectance ratio of the extreme infra- 
red band (MSS 7) to the chlorophyll band (MSS 5) increases with plant 
density. This ratio has been reported as a good data transformation 
for assessing changes in scene vegetative biomass (Maxwell, 1975), 
which is closely related to LAI, The characteristic increase in 
this ratio throughout the lower LAI's, followed by a relatively flat 
response of the higher indices agrees with the general results of 
these empirical studies. 

6.3 Simulation for Satellite Radiance 

The induced atmospheric effects and conversion from canopy re- 
flectance to predicted satellite radiance values was achieved by 
executing the Turner model with this nominal data set for each can- 
opy reflectance data set. In the case of canopy model generated re- 
flectances, each row of tne Simulated Canopy Reflectance table 
constituted an input. Each individual field measurement acted as 
input for the empirical case. These basic units were used, rather 
than average reflectance values, in order to facilitate the calcula- 
tion of a covariance matrix in terms of radiance units. The mean 
reflectance values depicted in Table 2 could have been used in 
establishing the mean radiance values, however, this simple ap- 
proach ignores the variance associated with each reflectance data 
point. By translating each reflectance value into radiance units, 
the variance can be directly calculated. 

An additional step is needed in converting the field measure- 
ments into estimated signals at a satellite. As the field data was 
collected throughout a day, considerable sun angle effects are 
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FIGURE 53. RATIO OF BANDS 7/5 
AS A FUNCTION OF LAI. 
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ingrained in the data. This variation is not present in the model 
data because all of the model runs for each phonological stage were 
made for a single sun position. In order to make the empirical and 
model data sets compatable, a sun angle correction factor was applied 
to each empirical measurement which converted it into terms of the 
model's sun position. 

The procedure for sun angel correction involved the conversion 
of field measured reflectance into radiance units, while maintaining 
continuity of sun angle. The result of this operation is an esti- 
mated sensor signal keyed to the particular sun angle occuring at 
the time of measurement. This signal can then be translated into 
terms of another sun angle by applying a linear correction algorithm 
(Smith, Berry and Heimes, 1975). For example, the March model pre- 
dictions were made for a sun angle of 51 degrees zenith. A field 
data set taken at 1730 hours denotes a 35 degree zenith sun angle. 

The correction procedure would involve the execution of the Turner 
model using the measured reflectance values and specifying a 35° 
degree sun angle. This prediction could then be expressed in terms 
of the model's sun position by solving the following equation: 

N 51 O = (»5i"_350 * N 35 O) + S 51 O. 35 O (5) 


where, 

NgjO = corrected radiance for 51° sun angle 
NsgO = radiance for 35° sun angle 
“51°. 35° " niultiplicative correction factor for adjusting 
35° to 51° sun angle 
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510 35 


0 = additive correction factor for adjusting 35“ to 
51° sun angle. 


This correction procedure was applied to each field measurement, 
and the mean vector and covariance matrix were calculated for each 
plot. Table 3 and Figures 54 and 55 present the radiance data for 
both the model and empirical reflectance data sets. The interpreta- 
tion of the general trends in the data is similar to chose presented 
for the reflectance data in the previous section. The only differ- 
ences are subtle changes in the plot values due to the correction 
for varying sun angles. 


6.4 Analysis of Model Data 

In general, the ability of the canopy reflectance model to track 
the empirical measurements is good. In addition to the comparisons 
outlined in this section, a diurnal reflectance comparison between 
model predictions and actual fit measurements was made and des- 
cribed in a report by Smith, Berry and lleimes (1975). The results 
of this evaluation v/ere also positive. However, three aspects of 
the data warrant detailed discussion. 

The first aspect is the relative inaccuracy of the March model 
predictions for MSS bands 4 and 5. The source of error in these 
predictions is most likely a result of the model's inability to 
adequately deal with the pronounced rowing effect at low LAI's, and 
the strong contribution of a highly variable soil reflectance. This 
problem is addressed in the canopy model in an input parameter de- 
noting spatial dispersion of foliage elements. However, this para- 
meter was not determinec in this study. 
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Table 3. MODEL/FIELD SIMULATED RADIAfJCE AS A FUNCTION OF LAI. 


LAI MSS 4 MSS 5 MSS 6 MSS 7 MSS 7/5 


MARCH (MODEL) 

.5 2.700 

1.0 2.951 

1.5 2.805 

2.0 2.981 

2.5 2.874 

3.0 2.973 

3.5 2.987 

4.0 2.989 

4.5 3.055 

5.0 3.055 

APRIL (MODEL) 

3.5 3.456 

4.0 3.469 

4.5 3.504 

5.0 3.420 

5.5 3.479 

6.0 3.460 

6.5 3.452 

7.0 3.422 

7.5 3.421 

8.0 3.421 

MARCH (EMPIRICAL) 

1.31 3.467 

2.07 3.520 

4.06 3.655 


5.13 3.333 
5.36 3.218 
6.15 3.248 


1.922 

3.670 

3.217 

1.674 

1.922 

5.642 

5.146 

2.677 

1.714 

5.528 

5.056 

2.950 

1.840 

6.204 

5.683 

3.089 

1.684 

6.148 

5.649 

3.355 

1.746 

6.527 

5.983 

3.427 

1.756 

6.595 

6.050 

3.445 

1.738 

6.666 

6.116 

3.519 

1.781 

6.929 

6.356 

3.569 

1.781 

6.934 

6.364 

3.573 

1.986 

7.598 

7.051 

3.550 

1.995 

7.6711 

7.151 

3.584 

2.018 

7.808 

7.268 

5.602 

1.966 

7.515 

7.020 

3.571 

2.001 

7.705 

7.157 

3.577 

1.989 

7.633 

7.091 

3.565 

1.983 

7.588 

7.045 

3.553 

1.963 

7.446 

6.899 

3.515 

1.962 

7.442 

6.895 

3.514 

1.962 

7.422 

6.895 

3.514 

2.626 

6.888 

6.922 

2.636 

2.750 

5.574 

5.142 

1.870 

2.824 

6.057 

5.434 

1.924 

1.891 

7.286 

7.422 

3.925 

1.733 

7.776 

7.753 

4.373 

1.780 

7.898 

7.803 

4.384 


APRIL (EMPIRICAL) 














The second aspect is the pronounced descrepency In March between 
model and field determined reflectance as a function of leaf area 
index. This disparity is particularly apparent in the infrared bands. 
The empirical data indicates a declining trend in reflectance as LAI 
is increased, whereas the model trend is the opposite. Previous 
studies in this area (Johnson, 1975; Pearson and Miller, 1973) tend 
to support the model derived trend. An atypical field reflectance 
measurement could result through the highly variable soil reflectance 
contribution and differing percent ground covers circumscribed by 
the field of view of the LANDSAT field radiometer. 

The f i nal aspect i s the degeneracy of some of the model ' s de- 
rived covariance matrices. The covariances betv/een all bands v^ere 
unusually high, which resulted in unfeasible divergence calculations. 
The source of this effect is most likely ingrained in the model's 
treatment of the leaf and soil radiometric parameters. The general 
form of the SRVC model allows all of the scene compon“snt tc-be sto- 
chastic. However, for its application in this study, leaf transmis- 
sion and reflectance and soil reflectance were represented by constant 
parameters. This covariance between spectral bands was therefore 
highly correlated. 


7.0 CLASSIFICATION PROCEDURE 


The previous sections of this paper have described and developed 
the data base necessary to infer leaf area index from plant canopy 
spectral measurements. A field measurements program v/as designed to 
establish physical parameters used to calibrate and execute a canopy 
reflectance model. This model was then used to derive a data base 
expressing the direct relationship between LAI and canopy reflec- 
tance. A second model was employed to simulate atmospheric effects 
and predict the resultant relationship between LAI and radiance re- 
ceived at a satellite sensor. In a similar manner, the field meas- 
ured canopy reflectance for several evaluation plots was transcribed 
into predicted radiance values. This section presents a description 
of the technique used in classifying plots of unknown leaf area 
index. 

7.1 The Swain-Fu Distance Measure 

The data classification procedure used in this research is a 
hybrid of the two general approaches used in pattern recognition. It 
is similar to supervised classification in that a priori knowledge 
is assumed. The descriptive statistics for the various classes of 
LAI, however, were developed through computer modeling, rather than 
empirical derivation from ground truth training sets. The actual 
classification algorithm used, employs distance measures in the same 
manner as unsupervised clustering techniques. The classification 
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procedure involves comparing the spectral response of an area with 
unknown LAI to those responses of known LAI. The divergence be> 
tween each unknown model response pair is calculated, and the un- 
known point is assigned an implied LAI equal to that of the "closest" 
model response. Throughout this discussion "effects'* will refer to 
the model generated responses used to classify the unknown LAI re- 
sponses of "plots". 

The most familiar distance or divergence measure is Euclidean 
distance. This technique employs only the mean vector to assess 
the point to point distance between two responses. It principally 
involves calculating a directionally independent difference between 
the means of two responses. In mathematical terms, the Euclidean 
distance for two n-dimensional points (2(, YJ is: 


^Euclidean * ( j 


( 6 ) 


A variation of Euclidean distance considers the dispersion of 
data as well as the mean response vector in comparing the separation 
of points. This dispersion is approximated in terms of the "ellip- 
soid of concentration", which is calculated from the mean vector and 


covariance matrix. This geometric configuration can be conceptual- 
ized as replacing a "cloud" of individual data points* plotted in 
three-space, with a football centered about the mean vector* The 
surface of the football describes the bounds of the ellipsoid of 
concentration. In contrast. Euclidean distance considers only the 
epicenters of the footballs, while the more advanced measures judge 
the shape and orientation as well. 
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Figure 56 is a two-space portrayal of the Swain-Fu distance meas- 
ure developed at Purdue University (Swain » 1973). In terms of the 
distance shown Swain-Fu is given by: 

^Swain-Fu Dj^ + Dg (7) 

In narrative terms, the divergence is calculated by weighting the 
sample Euclidean distance (O^g) by the inverse of the spread of the 
data along the axis connecting the centers of th*; clusters (Dj and 
Dg). In terms of th mean vectors and covariance matrices associated 
with the clusters, tha distance is expressed as, 


Swain-Fu 


(Ck.,)’^*(Ci.k)‘^ 



where , 


C, = tr 
1 ,k 

(Hi - UkHMi 

■V’l 

(9) 

^k.l = **' 



(10) 

tr(N) = trace of matrix N 




E * covariance matrix for cluster 
'-n 

U = mean vector for cluster n 

-fi ■ 

7.2 Covariance Matrix Modification 

the original intention of this process was to use the actual 
covariance matrices derived from the field met jurement for the 
"plots” , and those developed by the stochastic canopy model for 
the "effects". However, several of the model generated covariance 



matrices proved degenerate, and yielded undefined Swain-Fu distances 
An explanation as to the probable cause of the degeneracy is given 
in Section 6.4. 

The covariance matrices were modified by replacing the hyper- 
ellipsoid clusters with hyper-spheres of equal volume. Figure 57 
is a two-dimensional schematic of the important factors in this 
technique. The equation for the surface area of an ellipse is, 

A * irab (11) 

where, a and b = major and minor axes, respectively. The surface 
of a circle is calculated by, 

A = irr^ (12) 

where, 

r = radius. 

Principal Component theory shows that the axes of any elliptical 
form are the eigenvalues of the square matrix defining the form 
(Xj and Xg in the schematic). Thus, equation (11) can be rewrite 
ten as, 

A = irXj^Xg (13) 

Setting equation (11) and (12) equal to each other, and solving for 
r, results in, 

2 

irr = irXj^Xg 

r = (XjAj)'* (14) 

This final equation a^ the solution for the radius of a circle 
having the same surface area of a given ellipse, in terms of its 
eigenvalues. This same relationship, expressed in terms of 
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FIGURE 56. SWAIN-FU DISTANCE MEASURE. The 
Swain-Fu distance measure utilizes simple 
cluster distance (D 12 ) cluster dis- 

persion (Dj and Dg); 

(J) 




hyper-geometry, is. 


'^n-space " 


* \ * 


* X 1 

W 


1/n 


( 15 ) 


The derived hyper-sphere can then be substituted for the actual 
ellipsoid of concentration by replacing the actual covariance ma- 
trix with a diagonalized matrix based on the radius of the circle. 

This modified Sv/ain-Fu distance measure forms the foundation of 
the classification procedure used in this research. Program SUFI) 
(Appendix C) is the FORTRAN coded version of the computer algorithm 
developed. The program consists of four phases: 1) input and data 
modification: 2) weighted distance calculation; 3) classification; 
and <?) output. The first phase allows for the input of the scaled 
"effects", "plots" with unknown leaf area indices, and specific pro- 
gram parameters for control. Data modification can be performed 
on both the mean vectors and covariance matrices. The divergence 
calculation phase solves for the Swain-Fu distance measure for each 
"effect/plot" pair. This procedure first evaluates expressions (4) 
and (5) for the weighted one-way distances between clusters. These 
values, termed C- factor distances, are then substituted into equa- 
tion (3) to develop the overall divergence measures. The third 
phase of the routine determines the "effect/ plot" pair with minimal 
divergence. Leaf area index classification of the "plot" is achieved 
by its association with the closest "effect". The final phase out- 
puts the results in the forms of a C- factor matrix, Swain-Fu dis- 
tances, and a printer plot of divergence. Details of the specific 
procedures are contained in the SUFU program listing presented in 
Appendix C. 


8.0 CLASSIFICATION RESULTS 


The basic classification of the field plots by the SUFU rou- 
tine involved a 4-dimensional feature space, with the response in 
each of the LANDSAT bands forming a single feature vector. Clas- 
sification of the data was made by considering only Euclidean 
distances between the “effects" and the "plots" data clusters. Two 
additional approaches were also executed, in which the Swain- Fu 
distance was used. The first considered two raw data preprocessing 
transformations: 1) spherical elipsoid of concentration; and 2) cor- 
rection for sun angle effects. The final approach utilized the pre- 
processed data to form new feature vectors. These derived feature 
vectors include: 1) the ratio of bands 7 to 5; 2) the ratio of 
bands 6 to 4; and 3) the simultaneous consideration of the 7 to 5 
and 6 to 4 ratios. 

The field measured leaf area index for each of the field "plots" 
used in this study is identified in Table 4. During the March meas- 
urement period, “plot" 2 had an unusually high LAI of 4.60. "Plots" 

1 and 3 are more nearly normal, having LAI's of 2.07 and 1.31, re- 
spectively. The April measurements shov/ "plot" 3 as having a LAI of 
6.15, with "plots" 1 and 2 being nearly equal at 5.13 and 5.36, re- 
spectively. 
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Table 4. 

FIELD MEASURED LEAF AREA INDEX 




DATE 

PLOT 

LAI 

RANKING 

March 

i 

2.07 

2 


2 

4.06 

' 1 ,_, , 


3 

1.31 

3 

April 

1 

5.13 

;.-.3- . 


2 

5.36 

2- 


3 

6.15 

1 



3 
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8.1 Raw Data Classification 

In the raw data classification scheme, the SUFU program was used 
In a manner which approximated a Euclidean distance measurement of 
divergence. The mean vectors for the “plots" were developed by 
averaging the radiometer measurements taken throughout a day. This 
has the Impact of tacitly disregarding any sun-angle effects on the 
spectral signatures. The mean vectors for the "effects" consisted 
of the model predicted responses. 

In order to force the SUFU routine to act as a Euclidean dis- 
tance classifier, a common covariance matrix was used for all "plots" 
and "effects". The matrix was diagonalized, with the diagonal ele- 
ments set to unity. This operation has the effect of containing 
the dispersion of the data about each mean vector, v/hich removes any 
classification information based on differential data dispersion. By 
default, all of the classification intelligence became embedded in 
the mean vector positioning. 

Figure 58 shows the results of the raw data classification using 
the measured and modeled reflectance values. The March assignments 
show very strong classifications. This is indicated by the minimum 
normalized divergence for each case being less than .5. In des- 
criptive terms, a .5 normalized divergence denotes a cluster which 
Is half the distance between the "plot" and the farthest away "ef- 
fect". At this level, an appreciable difference between the closest 
"effect" and the other "effects" appears to exist, less pronounced 
classifications occur for the April reflectance data, as shown by 
the minimal normalized divergences being about .6. 


NORMALIZED DIVERGENCE 
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0.5 1.5 2.5 3.5 4.5 3.5 4.5 5.5 6.5 7.5 

LEAF AREA INDEX 

FIGURE 58. RAW REFLECTANCE DATA DIVERGENCE PLOT. 



The classification accuracy of the raw reflectance data appears 
best during the April period. The program predicts a LAI of 4.5 
for all three “plots", which is close to their 5.5 average LAI. The 
March predictions, however, are less accurate and more difficult to 
interpret. The assigned LAI for both "plots" 1 and 2 is 1.0, with 
"plot" 3 indicating an LAI of 5.0. The assignment for "plot" 1 is 
relatively near the measured 2.07 LAI. The classifying of "plots" 

2 and 3, on the other hand, are nearly the reverse of the measured 
4*60 and 1.31 LAI for the respective plots. 

Figure 59 identifies the LAI classification based on the unal- 
tered radiance data predicted by the atmospheric model. The classi- 
fication results and divergence graphs are nearly the same as those 
of the previous case. The notable exception occurs for "plot" 1 
during the April period. The classification of this "plot" changes 
slightly from 4.0 to 4.5 LAI. The general similarity of these clas- 
sifications Indicates that the relative positioning of the mean 
vectors -are somewhat insensitive to atmospheric effects. 

8.2 Preprocessed Data Classification 

Figure 1 of Appendix B portrays the results of the SUFU classifi- 
cation based on the covariance weighted distance measure. It is 
apparent that the consideration of covariance in the divergence cal- 
culation using surface reflectance data had no impact on LAI classi- 
fication, and only minimal effects on the normalized divergence 
plots. The result indicates that the covariance matrices for both 
the "plots" and "effects" tend to be similar, and without any d(Mn- 
inating lobes. Subsequent review of the covariance matrices of the 
data, support this contention. 
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A sharp contrast is noted for the classification based on 
radiance data, using the spherical approximation of the actual 
covariance matrices. Figure 3 of Appendix B summarizes these 
results. Whereas the Euclidean classifier tagged the March "plot" 

2 as having an LAI of 1.0, the covariance weighted procedure iden- 
tified it as having an LAI of 2.0. In a similar manner, April 
“plots" 1 and 3 are classified as having LAI's of 5.0 and 4.5, 
respectively, by Euclidean distance. The data dispersion weighted 
technique, however, identifies these same plots as both having an 
LAI of 4.0. 

These descrepencies are most likely the result of the incorpor- 
ation of spherical approximations of the data dispersion. Another 
possible explanation is that the dispersion of the data is selec- 
tively altered by atmospheric effects. As noted earlier, the mean 
vector positioning, however, is relatively unaltered. The within 
wavelength band variance has a possible physical explanation in the 
effects of selective absorption and scattering affecting the path 
radiance component of the satellite signal. 

In addition to the data preprocessing for spherical covariances, 
a case was executed in which the "plot" data for both periods were 
corrected for sun angle effects. A linear correction algorithm, 
specifically developed for this data, was used. The procedure first 
involved relating the time of day for each canopy reflectance meas- 
urement to the appropriate sun angle. Once the sun angle for each 
field measurement was detemined, the corresponding mean radiance 
prediction could be corrected to the base sun angle used in the model 
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derivation of the "effect". The theoretical foundation and details 

t I 

■ } 

of this correction technique are presented in a report by Smith, 

Berry and Heimes (1975). 

The classification results, shown in Figure 3 of Appendix B are 
nearly the same as those for the uncorrected data. The divergence 
plots are also very similar. This outcome appears to demonstrate 
that sun angle effects, over the range of sun angels involved in 
this study, have a negligible effect on LAI classification. 

8.3 Derived Feature Vector Classification 

To this point, the preprocessing of data has concentrated on 
meaningful transformations of the responses themselves. Another 
technique frequently used in remote sensing is data transformations 
designed to derive new feature vectors. These include ratioing , 

4 - 

bands, adding or subtracting them, and applying arithmetic func- 
tions to stretch and distort the original data. The purpose of 
these operations is to identify new feature spaces which enhance 
particular characteristics of a scene. For example, a relatively 
low response in band 4, when compared to band 7, might provide the 
same amount of information as the four band signatures in identi- 
fying vegetated areas. 

Figures 4 through 6 of Appendix B report the results of using 
band ratios of 7/5, 6/4, and 7/5 with 6/4, The data used in these 
feature vector transformations utilized spherical covariance ma- 
trices and sun angle corrected mean vectors. The selection of these 
transformations was based on the work by 6. Johnson (1976) in asses- 
sing plant biomass from LANDSAT data. The results, for the most 
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part appear disappointing. Numerous departures from the previous 
classifications are noted, and the normalized divergence plots of- 
ten display minimal deflections. This erratic behavior is most 
likely caused by the severe averaging effect of the ratioing and 
the minimal number of feature vectors for classification. Of the 
three transformations, the 6/4 ratio is the closest to the original 
classifications, but underestimates the March leaf area indices. 

8.4 Sumnary of Classification Results 

A suninary of the classification results for the radiance data of 
this study is reported in Table 5. The table is constructed in ac- 
cordance with the major classification categories. The row entries 
are organized according to phenology stage and field plot number. 

The column organization reflects the three classification techniques 
1) raw data using Euclidean distance; 2) corrected data using Swain- 
Fu divergence; and 3) ratioed data, also using the Swain-Fu measure. 

The best overall results were achieved by the simple classifica- 
tion of mean spectral responses using a Euclidean distance classi- 
fier. The introduction of data dispersion information in 
classification, through the consideration of spherical covariance, 
actually degraded the results. Preprocessing the data for sun angle 
only slightly improved the classification accuracy. The use of 
derived feature vectors resulted in the worst classifications, and 
consistently understated the leaf area indices for the March period. 

The ability of the procedure to differentiate between conditions 
of average high plant density and low density is apparent. During 
the low density tillering stage in March the average deviation of 



B 
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Table b. SUMMARY OF CLASSIFICATION RESULTS. 
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predictions based on technique I v«s .15. The worst technique 
(ratio of bands 7/5) yielded an average deviation of 1.81 for this 
same period. The average deviations associated with the more dense 
heeding stage in April proved to be more consistent between tech- 
niques, resulting in average deviations from .22 to 1.38. 


9.0 CONCLUSION 


The fundamental purpose of this research was to investigate the 
feasibility of using abstract computer modeling to infer inherent 
scene characteristics from remote sensing data. The principle hy- 
pothesis of the study was: 

- Plant canopy reflectance and atmospheric modeling 
can be used to infer intrinsic scene geometry var- 
iables from spectral measurements. 

The specific geometric variable under study was the leaf area 
index (LAI) of a wheat canopy. A stochastic plant canopy reflec- 
tance model was used to develop a simulated reflectance data base 
which tracked the induced effects of changes in the LAI variable. 

An atmospheric model was employed to translate the surface reflec- 
tance predictions into simulated satellite signals. The results of 
these operations were two data sets responding to systematic chan- 
ges in LAI: one in terms of scene reflectance » and the other iden- 
tifying radiance. Classification of measured field reflectance and 
corresponding simulated radiance was achieved by determining the 
minimum divergence between the model spectral signatures of known 
LAI* and the signatures of the responses to be identified. The 
leaf area index for an unknown area was assigned the same value as 
that of its most similar model derived signature. 

In partial support of this approach* an intensive field meas- 
uranent program was conducted. These data were utilized in the 


study in two ways. First* the data acted as input to the canopy re- 
flectance model and allowed for the calibration of the model's per- 
formance. Tlie second major use was in the evaluation of the hypotheses 
under investigation. 

9.1 Hypotheses Results 

In addressing the general hypothesis of using modeling to infer 
scene geometry from spectral data, two sub-hypotheses were madet 

- Wheat canopy reflectance can be predicted by computer 
modeling, both as a function of low and high plant 
densities, and as a function of sun angle. 

- Low wheat density and high wheat density categories 
can be inferred from model data sets. 

The results of the investigation into the first sub-hypothesis 
demonstrates that reflectance modeling is applicable to wheat can- 
opies. Section 6.2 describes this investigation. Figure 50 and 
Tables 2 and 3 summarize the results of this study. It was found 
that the general agreement between model and field spectral signa- 
tures for both the March (low LAI) and the April (high LAI) periods, 
were excellent. The March period was less accurate, most likely as 
a result of the pronounced rowing effect of the canopy, and an ini- 
tial assumption made by the author which constrained the model's 
operation. The ability of the reflectance model to track sun angle 
effects was also apparent. A detailed account of this result is re- 
ported in a paper by Smith and this author (1975). 

The results of the second sub-hypothesis is reported in Section 
8.0 and Table 5. During the low density stage in March the average 
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measured LAI was 2.48 while the average model prediction* using the 
best techn1^ue, Inferred a LAI of 2.33. A comparison of April LAI's, 
noted an actual average of 5.55, and an Inferred Index of 4.67, 
using the same classification technique. 

The general hypothesis Is Inconclusive when applied to Individual 
plots. This results primarily from the limited number of field plots 
measured In this study as discussed In Section 8,4. 

9.2 Other Major Contributions of the Study 

In addition to the major efforts of applying and validating a 
reflectance model to wheat canopies, and the development of a meth- 
odology for Inferring leaf area index, several other contributions 
of this study should be noted. An extensive field measured data 
base was collected which Incorporates the primary environmental and -| 

intrinsic scene variables, with actual canopy reflectance, for the 
four major phenology stages of wheat. This data base Is further 
extended by the addition of model derived canopy reflectance as 
functions of both sun angle and leaf area Index. 

The empirical study was also valuable In that several unique 
data collection techniques were developed. In particular, were the 
design and construction of an instrument for field measurement of 
individual leaf transmission (Section 3.2.1), and two rapid, in situ. 
field techniques for assessing leaf angle distributions (LAD) Sec- 
tions 3.2.2, 3.3, 4.0 and 5.0. 

Two major mathematical contributions were made. The sensitivity 
analysis of the Fredholm technique (Section 4.1) describes a proce- 
dure for evaluating the magnitude and direction of the sensitivity ( ) 


of an individual input vector on an indivisible output vector. The 
Fourier technique required the development of a mathematical proce- 
dure for convoluting two orthogonal distributions of angles, in 
order to estimate the true 3-space leaf angle distribution of a 
plant canopy (Section 4.2,2). 

The final contribution of the study involves an evaluation of 
several data transformations and classification algorithms, as to 
their effect in classifying leaf area index from spectral data 
(Section 8.0). The specific procedures investigated were; 1} 
Euclidean and covariance weighted classifiers; 2) sun angle correc- 
tions; and 3) derived feature vectors. In conjunction with the 
first procedure a technique was implemented which transforms the 
ellipsoid of concentration of a data cluster into a hyper-sphere 
approximating the data dispersion information. 

9.3 Future Research 

The most apparent area of future research associated with this 
studiy, is the evaluation of the procedure under more variable con- 
ditions and with actual aircraft or satellite data. In a study of 
this type, an operational test could be designed in which extensive 
ground measurements of LAI vould serve as evaluation points. The 
hypotheses of such a study might include the principle hypothesis 
outlined in this report and a sub-hypothesis concerned with the 
extent of varietal or physical condition variations. 

A second area of potential study is ingrained in the refinement 
and extension of the orocedure. This research might be concerned 
with the modeling approaches and assumptions. In addition, a more 
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rigorous investigation into classification techniques is warranted. 
Finally* the basic technique could be employed to determine the 
feasibility of inferring other intrinsic scene parameters, e.g., 
individual leaf reflectance and transmission. 
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APPENDIX A 


Radiometric and Geometric Field Data 

The principle field data collected by TAMU/CSU for the canopy 
modeling effort consists of periodic canopy reflectance, intensive 
leaf area index (LAI) measures, extensive LAI estimates, individual 
leaf transmission measurements, and canopy geometry photos. Rela- 
tively complete data sets are available for March 20, 1975 (Til- 
lering Stage, TAMU), April 23, 1975 (Jointing Stage, TAMU/CSU), and 
May 20, 1975 (Heading Stage, TAfiU). Less complete data sets were 
collected on November 24, 1974 (Winter Tillering Stage, TAMU/CSU) 
and June 26, 1975 (Ripening Stage, TAMU/CSU). The following table 
summarizes the data set. A more detailed presentation of the March 
and April data sets is included in the remaining parts of this 
section. 
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Table A-1. FINNEY COUNTY DATA SUMMARY, The field measurements of 
the radiometric and geometric parameters of the canopy 
were made for four phonological stages. Diurnal can- 
opy reflectance was collected for each stage. The 
diffuse to direct i radiance ratio and soil reflectance 
were sampled periodically for each data. 


March 20, 1975 

Tillering Stage 

Field 416 

-Canopy Reflectance: 





Plot 1 

Plot 2 

Plot 3 

Time: 

1100 hrs. 

WT 


1030 


1145 

1130 


1115 


1300 

1245 


1230 


1400 

1345 


1330 

-Leaf Area Index: 

2.07 

4.06 


1.31 

-Canopy Gecmietry: Fredholm Field Photos 



-Leaf Transmission: Not Taken 




-10" LAI Plots: Field 

367 369 

370 

414 

421 

Plot 1 

2.22 2.15 

1.48 

8.53 

4.54 

Plot 2 

2.78 2.43 

8.45 

9.17 

3.60 


April 23, 1975 

Jointing Stage 

Field 416 

-Canopy Reflectance: 

Plot 1 

Plot 

_2 

Plot 3 

Time: 

1000 hrs. 

TMT 


1115 


1130 

1145 


1200 


1315 

1345 


1400 


1715 

1730 


1800 

-Leaf Area Index: 

5.13 

5.36 


6.15 

-Canopy Geometry: Fourier Field Photos 



-Leaf Transmission: Not Taken 




-10" LAI Plots: Field 

367 369 

370 

414 

421 

TOTl 

2.22 2.15 

1.48 

8.53 

4.54 

pWT 

2.78 2.43 

8.45 

9.17 

3.60 



III. May 20, 1975 

Heading Stage 

Field 416 

-Canopy Reflectance: 

Plot 1 

Plot 

2 

Plot 3 

Time: 

ti5?rFrs. 

mr- 

1045' 


1100 

1115 


1130 


1200 

1215 


1245 


1300 

1315 


1345 

-Leaf Area Index; 

4.11 

5.32 


6.04 

-Canopy Geometry: Fourier Field Photos 



-Leaf Transmission; Green, Yellowing 

, Dead 



-10" LAI Plots: Field 

367 369 

370 

414 

421 

WoTl 

3.82 1.83 

3.12 

5.65 

1.80 

Plot 2 

3.22 5.68 

9.76 

7.64 

2,16 


IV. June 26, 1975 

Ripening Stage 

Field 416 

-Canopy Reflectance: 

Plot 1 

Plot 

1 

Plot 3 

Time: 

-- 

1115 


1000 


-- 

1200 


1115 


-- 

1245 


-- 


— 

1300 


— 

-Leaf Area Index; 

1.79 

2.17 


2.04 

-Canopy Geometry: Fourier Field Photos 



-Leaf Transmission: Dead 





-10"LAI Plots: Field 

367 369 

370 

414 

421 

FT^l 

.78 1.15 

.94 

1.79 

1.02 

Plot 2 

.85 2.00 

3.14 

2.63 

1.08 
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Table A-2. MARCH LEAF AREA INDEX AND DESCRIPTIVE PARAMETERS. 


March 20, 1975 Tillering Stage Field 416 


Santana Wheat 


Crop Type 
Height 
Chlorotic 
Weeds 

Soil Condition 
Wind 

Leaf Area Index 


8-9 cm. 

Green Foliage 
0 % 

Moist 

12-15 mph NW 


(10" drill; EW) 


PLOT 1 
2.07 

58.10 gm 

803.00 

0.00 

803.00 

8.64 

0.00 

8.64 

39.39 cm^ 
0.00 
0.00 

37.39 


PLOT 2 PLOT 3 
4.60 1.31 

39.90 

1372.00 

0.00 

1372.00 

11.00 

0.00 

11.00 

60.09 

0.00 

0.00 

60.09 


Dry Weight (2* X 2' Plot) 

Number of Tillers (2‘ X 2* Plot) 
Live 
Dead 
Total 

Average Tillers/Plant 
Live 
Dead 
Total 

Average Lead Area/ PI ant 
Green 
Yellow 
Dead 
Total 


PROBABILITY 
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FIGURE A-1. LEAF ANGLE DISTRIBUTION FOR MARCH 



Table A-3. MARCH RADIOMETRIC DATA 


DATE 

TIME 

CROP AND LOCATION 

- PLOT NUMBCR 

032075 

1052 

HHEAT KS. 

416-1 

032075 

1055 

WHEAT KS. 

416-1 

032075 

1057 

WHEAT R5. 

416-1 

032U75 

I05B 

WHEAT K5. 

416-1 

032UTS 

1136 

HHEAT KS. 

416-1 

032075 

1135 

WHEAT KS. 

416-1 

032075 

IMO 

wheat KS. 

416-1 

032075 

1142 

WHEAT KS. 

416-1 

032075 

0104 

wheat KS. 

416-1 

032075 

0105 

wheat KS. 

416-1 

012075 

0104 

wheat KS, 

416-1 

012075 

0107 

WHEAT KS. 

416-1 

0320 75 

0155 

wheat KS. 

416-1 

032075 

0156 

WHEAT KS. 

416-1 

032075 

0155 

wheat KS, 

416-1 

032075 

0159 

wheat KS, 

416-1 

032075 

04.15 

wheat ns. 

416-1 

032075 

0*36 

WHEAT KS. 

416-1 

032075 

0437 

wheat KS. 

416-1 

032075 

0441 

wheat KS. 

416-1 

037075 

1044 

wheat KS, 

416-2 

032075 

1045 

WHEAT KS. 

416-2 

032075 

1046 

wheat KS. 

416-2 

032075 

1047 

WHEAT KS, 

416-2 

032075 

1127 

wheat KS, 

416-2 

032075 

1127 

wheat KS. 

416-2 

032075 

1125 

Wheat ks. 

416-2 

0320 75 

1179 

wheat KS. 

416-2 

032075 

1252 

wheat ks. 

416-2 

032075 

1254 

wheat ks. 

416-2 

032075 

1256 

wheat ks. 

416-2 

032075 

* . ^9 > 

wheat ks. 

416-2 

032075 

r •? 

wheat ks. 

416-2 

0320 75 

01*2 

WHEAT KS. 

416-2 

032075 

01*3 

wheat ks. 

416-2 

032075 

0144 

wheat ks. 

416-2 

032075 

0115 

WHEAT KS. 

4I6-? 

032075 

5335 

wheat KS. 

416-2 

032075 

0339 

WHEAT KS. 

416-2 

032075 

0 341 

WHEAT KS. 

416-2 

032075 

ili3« 

wheat ks. 

*16-3 

032075 

1035 

WHEAT KS. 

416-.1 

037075 

1036 

wheat KS. 

416-T 

032075 

1037 

wheat Kb. 

416-3 


orientation 

REFLECTANCCk HANOI 

HANn2 

R4N03 

-RANP* 

OFF RON 

.074 

.081 

.742 

.335 

ON ROW 

.075 

.079 

.?5* 

.320 

ON ROW 

.977 

.075 

,72* 

.793 

OFF ROW 

.074 

.07* 

.719 

.310 

OFF ROW 

.06.3 

.070 

.717 

.315 

ON ROW 

.069 

.075 

.722 

,337 

ON ROW 

.0*3 

.069 

.727 

.301 

OFF RON 

.066 

.07? 

.708 

,307 

OFF RON 

.065 

.060 

.715 

.294 

ON ROW 

.060 

.060 

.190 

,250 

On row 

-,060 

.050. 

,21* 

.31* 

off ROW 

.063 

.065 

.715 

.304 

ON ROW 

.05? 

.070 

.705 

.29* 

off now 

.076 

.075 

.716 

.30T 

OFF ROW 

,070 

.074 

.711 

.291 

ON ROW 

.070 

.061 

.’05 

.30* 

off row 

.074 

.075 

.*42 .. 

-.311 

ON row 

.061 

.077 

.710 

.32* 

on row 

.059 

.073 

.744 

.3*3 

off row 

.07* 

.075 

.262 

,3*T 

ON ROW 

.050 

.066 

.244 

.333 

OFF ROW 

.052 

.064 

.759 

.359 

OFF ROW 

.079 

.077 

.776^ 

,355 

ON BOW 

.075 

.074 

.759 

,310 

OFF ROW 

.051 

.069 

.750 

.335 

ON ROW 

.060 

.071 

.239 

.314 

ON ROW 

.0*8 

.069 

.235 

.33* 

OFF ROW 

.050 

.070 

.250 

.335 

OFF ROW 

.065 

.067 

.74* 

.330 

ON ROW 

.067 

.075 

.735 

.333 

ON ROW 

.065 

.06* 

.725 

.303 

OFF ROW 

.065 

.070 

.241 

.337 

ON ROW 

.063 

.060 

.239 

.333 

OFF ROW 

.065 

,07l 

.239 

.32“ 

ON ROW 

.065 

.074 

.222 

*314 

OFF ROW 

.069 

.071 

.725 

.311 

ON ROW 

,050 

.09? 

.753 

,360 

off row 

.073 

.ns? 

.74* 

.333 

On row 

.071 

.05* 

.729 

.327 

ofk row 

.071 

.0«? 


.32*. 

On row 

.070 

.077 

.300 

,510 

ON ROW 

.07* 

.050 

.'21 

,*7o 

On how 

.070 

.070 

.’70 

,«50 

OFF ROM 

.070 

.OHO 

.300 

.**» 



ffic<?rsy5r>‘: 


Table A-3- cent. 


032075 

1114 

WHEAT 

KS. 

416-3 

032075 

1115 

WHEAT 

KS. 

416-3 

032075 

1116 

WHEAT 

KS. 

4X6-3 

032075 

1117 

WHEAT 

K5, 

416-3 

032075 

1236 

WHEAT 

KS, 

416-3 

032075 

1237 

WHEAT 

KS. 

416-3 

032075 

1230 

wheat 

KS. 

416-3 

032075 

1239 

wheat 

KS. 

416-3 

032075 

0126 

wheat 

KS. 

416-3 

032075 

0128 

WHEAT 

KS. 

416-3 

032075 

0130 

WHEAT 

KS. 

416-3 

032075 

0131 

WHEAT 

KS, 

416-3 

032075 

0423 

wheat 

KS. 

416-3 

032075 

0424 

WHEAT 

KS. 

416-3 

032075 

0425 

wheat 

KS. 

416-3 

032075 

0426 

WHEAT 

KS. 

416-3 


OFF ROW 

.070 

.070 

.300 

.43" 


ON ROW 

,070 

.070 

.300 

,450 


ON ROW 

. .060 

.080 

- .290 

.460 


OFF ROW 

.080 

.080 

.310 

.430 


ON ROW 

.060 

.050 

.240 

,380 


OFF ROW 

.060 

.050 

.260 

,380 


OFF ROW 

.050 

.050 

.260 

.380 


OFF ROW 

.050 

.050 

.240 

.390 


OFF RON 

.060 

.060 

.270 

.,390 


ON ROW 

.060 

.050 

.250 

.400 


ON ROW 

.060 

.050 

.250 

.400 


OFF PON 

.060 

•«60 

_.?60 

.3»0 


OFF ROW 

.073 

.080 

.300 

,440 


ON ROW 

.070 

.070 

.*70 

,4?0 


ON RON 

.070 

.070 

.290 

.420 


OFF ROW 

.080 

.090 

.290 

.440 
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Table A-4. APRIL LEAF AREA INDEX AND DESCRIPTIVE PARAMETERS. 



April 23, 1975 

Jointing Stage 

Field 416 


Crop Type 
Heigh c 
Chlorotic 
Weeds 

Soil Condition 
Wind 


Satanta Wheat (10" drill; EW) 
28-35 cm 
2-5% Yellowing 
0 % 


Dry 

Calm 



PLOT 1 

PLOT 2 

PLOT 3 

Lear Area Index 

5.13 

5.36 

6.15 

Dry Weight (2‘ X 2' Plot) 

142.24 gm 

148.70 

207.63 

Number of Tillers (2‘ X 2* Plot) 

Live 

907.00 

1179.00 

931.00 

Dead 

52.00 

75.00 

33.00 

Total 

959.00 

1254.00 

964.00 

Average Tillers/Plant 

Live 

8.80 

10.80 

8.40 

Dead 

.60 

1.60 

.80 

Total 

9.40 

12.40 

9.20 

Average Leaf Area/Plant 

Green 

100.84 cm^ 

83.48 

89.70 

Yellow 

33.21 

42.50 

19.16 

Dead 

40.99 

32.02 

31.49 

Total 

175.04 

158.00 

140.35 



FSOBABILm 
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Table A-5. APRIL RADIOMETRIC DATA. 



DATE 

time 

CROP AND LOCATiDN 

PLOT NUMBER 


042375 

100<» 

MHEaT ks. 

416-1 


042375 

1013 

wheat ks. 

416-1 


042375 

1015 

whe«t ks. 

416-1 


042375 

1016 

wheat ks. 

416-1 


042375 

1120 

WHEAT KS, 

416-1 


042375 

1130 

WHEAT KS. 

416-1 

042375 

1133 

wheat ks. 

416-1 

a 

042375 

1135 

WHEAT KS. 

416-1 

^ o 

042375 

0120 

WHEAT KS. 

416-1 


042375 

0124 

wheat KS. 

416-1 


042375 

0125 

WHEAT KS. 

416-1 


042379 

0127 

WHEAT KS. 

416-1 

50 

042375 

0521 

WHEAT KS, 

416-1 

0 

042375 

0523 

WHEAT KS. 

416-1 


042375 

0525 

wheat ks. 

416-1 

0 

042375 

0527 

wheat ks. 

416-1 

'r S 

042375 

1042 

WHEAT KS, 

416-2 


042375 

1044 

WHEAT KS. 

416-2 


042375 

1046 

WHEAT KS. 

416-2 

042375 

Iu4h 

WHEAT KS. 

416-2 


042375 

1146 

WHEAT KS. 

416-2 


042375 

1I4B 

wheat ks. 

416-2 


042375 

1150 

WHEAT KS. 

416-2 


042375 

1153 

WHEAT KS. 

416-2 


042375 

0140 

WHEAT KS. 

416-2 


042375 

0143 

WHEAT KS. 

416-2 


042375 

0146 

wheat ks. 

416-2 


042375 

0140 

wheat ks. 

416-2 


042375 

0538 

WHEAT KS. 

416-2 


042375 

0540 

WHEAT KS, 

416-2 


042375 

0541 

WHEAT KS. 

416-2 


042375 

0543 

wheat ks. 

416-2 



ORIENTATION 

REFLECTANCE* RANDl 

BAN02. BANDS BAN04 

ON 

ROW 

.037 

.023 

.?B4 

.407 

ON 

ROW 

.036 

.019 

.251 

.365 

OFF 

ROW 

.050 

.036 

.252 

.393 

OFF 

ROW 

.041 

.026 

.762 

.357 

, ON 

ROW 

. .043. 

. .027 

.280 . 

.434 

ON 

ROW 

,047 

*030 

.»T2 

,437 

OFF 

ROW 

.050 

.035 

.254 

.394 

OFF 

ROW 

.049 

..03B. 

.,?24 

.369 

OFF 

ROW 

.070 

.051 

.240 

.346 

OFF 

ROW 

.055 

.043 

,222 

.367 

ON 

ROW 

.053 

-.036 

-.292 . 

.462 

ON 

ROW 

.052 

.036 

.»77 

.443 

OFF 

ROW 

.053 

.024 

.279 

,395 

off 

ROW 

.05ft 

.020 

.»»B 

.404 

On 

ROW 

.051 

.033 

.753 

.389 

off 

RON 

.05? 

.031 

,270 

.475 

On 

ROW 

.039 

.020 

.284 . 

.450 

On 

ROW 

.039 

.025 

.*83 

.434 

off 

ROW 

,03R 

.023 

.243 

.350 

OFF 

ROW 

.041 

.026 

.254 

.371 

ON 

ROW 

.043 

.031 

.*77 

.413 

ON 

ROW 

.042 

.029 

.796 

.423 

OFF 

ROW 

,041 

.029 

.*21 

.340 

OFF 

ROW 

.044 

.031 

.247 

.346 

ON 

RON 

,047 

.035 

.300 

.437 

ON 

ROW 

.051 

.037 

.*94 

.441 

OFF 

ROW 

,049 

.03S 

.*56 

.356 

OFF 

ROW 

.047 

.032 

.*51 

,372 

ON 

ROW 

.053 

, 02 s 

.343 

,577 

ON 

ROW 

.057 

.029 

.352 

.571 

OFF 

ROW 

.04R 

.024 

.294 

.429 

OFF 

ROW 

.056 

.021 

.*93 

.433 


Table A* 5- cent 


0423 Tb 

Hill 

042375 

1113 

042375 

1114 

042375 

Ills 

0423/5 

1203 

0423/5 

1207 

042375 

1210 

042375 

1212 

042375 

0207 

042375 

0209 

042375 

0211 

042375 

0212 

0423/5 

0554 

0423/5 

0555 

042375 

0557 

042375 

0559 


KHEAT hS. 
•(MEAT KS. 
WHEAT KS. 
WHEAT KS. 
WHEAT KS, 
WHEAT KS. 
wheat KS. 
wheat KS. 
WHEAT KS. 
wheat KS. 
WHEAT KS. 
WHEAT KS, 
wheat KS. 
wheat KS. 
WHEAT KS. 
WHEAT KS. 


416-3 
416-3 
. .416-3 
416-3 
416-3 
4J6-3 
416-3 
416-3 
416-3 
416-3 
416-3 
416-3 
416-3 
416-3 
. 416-3 
416-3 


ON HOW 
ON ROW 
- ... OFF ROW 
OFF ROW 
ON ROW 
ON ROW 
OFF ROW 
OFF ROW 
. ..ON ROW 
ON ROW 
OFF ROW 
OFF ROW 
ON ROW 
OFF ROW 
. _ OFF Row 
OFF ROW 


.047 ,034 

.033 .019 

• 036 . .033 

• 044 ,021) 

.04? ,02T 

•OSO ..034 
.045 ,03? 

.041 .024 

.050 .,036.. 

.046 .034 

.044 ,03T 

•OSO .034 
•050 .026 

•052 .02T 

•060. ,023 
^056 ,023 


.305 ,AS3 
.»30 ,347 

.21" .344 

."*3 .421 

.263 .40T 

.’91 .464 

,"80 ,393 

."30 ,3?3 

.30 1 ,450 

."53 .415 

."64 ,36k 

."86 ,40? 

•32T ,521 

.323 ,50? 

,3.37 ,47a 

•364 ,50A 


TO 

- 1 ^ 


-f ■ 
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Table A-6. AVERAGE DIURNAL CANOPY REFLECTANCE. Canopy reflectance 
is averaged over on-rov/ and off-rov/ set ups. 


March 20 

Band 4 

Band 5 

Band 6 

Band 7 

PLOT 1 
10!52 hrs. 

.076 

.078 

.235 

.315 

1136 

.065 

.072 

.217 

.316 

1304 

.062 

.059 

.209 

.299 

1355 

.067 

.070 

.210 

.299 

1635 

.067 

.075 

.240 

.329 

PLOT 2 
iW 

.075 

.071 

.260 

.230 

1127 

.074 

.070 

.244 

.333 

1252 

.066 

.070 

.237 

.326 

1342 

.067 

.069 

.231 

.322 

1515 

.C74 

.086 

.242 

.335 

PLOT 3 
1033 

.071 

.083 

.298 

.475 

1114 

.055 

.075 

.300 

.443 

1236 

.055 

.050 

.250 

.333 

1331 

.060 

.055 

.258 

.390 

1623 

.068 

.073 

.288 

.430 


April 23 

r d 4 

Band 5 

Band 6 

Band 7 

PLOT 1 
1009 hrs. 

.023 

.027 

.249 

.381 

1128 

.047 

.033 

.258 

.409 

1302 

.058 

.058 

.258 

.405 

1738 

.054 

.029 

.270 

.416 

PLOT 2 
1042 

.039 

.025 

.266 

.401 

1146 

.043 

.030 

.261 

.381 

1340 

.049 

.035 

.275 

.402 

1738 

.054 

.025 

.321 

.503 

PLOT 3 
1110 

.040 

.025 

.262 

.391 

1203 

.045 

.031 

.266 

.395 

1407 

.049 

.035 

.276 

.408 

1754 

.055 

.025 

.338 

.519 



APPENDIX B 


Normalized Divergence Plots 

The divergence between th? model derived spectral signatures 
and the signatures of areas with unknown LAI were calculated 
using several different techniques. Figures B-1 through B-3 
show the normalized divergence plots for transformed data. Fig- 
ures b-4 through B-6 depict the normalized divergence for derived 
feature vectors. 
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FIGURE B-1. NORMALIZED DIVERGENCE USING SPHERICAL COVARIANCE MATRICES 

{REFLECTANCE DATA). 

PRECEIDING PAGE BLANK NOT FILMEXk 
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FIGURE B-2. NORMALIZED DIVERGENCE USING SPHERICAL COVARIANCE MATRICES 

(RADIANCE DATA). 
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0.5 1.5 2.5 3.5 4.5 3.5 4.5 5.5 6.5 7.5 

LEAF AREA INDEX 

FIGURE B-3. NORMALIZED DIVERGENCE USING SUN ANGLE CORRECTED RADIANCE. 
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FIGURE B-4. NORMALIZED DIVERGENCE USING RATIO OF BANDS 6/4. 
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FIGURE B-5. NORMALIZED DIVERGENCE USING RATIO OF BANDS 6/4. 
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APPENDIX C 


Computer Programs 

Five FORTRAN programs are presented in this appendix. The 
first three (PROP, THETA and CONVOL) are used in converting the 
measured angular bias of Fourier diffraction patterns, generated 
from field photographs, into estimates of the 3-space leaf angle 
distribution (LAD) of a plant canopy. Program ORTHOG utilizes 
digitized silhouettes of individual plants in calculating the 
distribution of angles for individual plants, and the overall LAD. 
Program SUFU is the classification routine developed in this study. 
All other programs used in the study are available through the 


references. 


Page intentionally left blank 
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Exhibit C-1. PROGRAM PROP. 


Program Name: PROP 
Subroutines Required: GRAPH 

Narrative: 

PROP accepts as input desitometer readings which wedge sample 
the diffraction pattern of an orthogonal view of a plant canopy. A 
plot of the distribution of leaf slopes contained in the original 
image is then generated. 

Control Card Input: 

Card 1 

Column 1-4 (14) (NRUN) Number of runs 

Column 5-10 (16) (NDATE) Date in 6 INTEGERS 

Column 11-20 (F10.5) (DERCT) Threshold value 

Column 21-30 (F10.5) (BADJ) Base adjustment for 

aperture 

Column 31-40 (F10.5) (DTEST) Minimum divergence 

test value 

Column 41 (II) (MTRAIL) Test for end of data 

(other than 0 for 
end of data) 

Cards 2 & 3 
Col umn 1-80; 

Card 2 (16F5.1) DATAE) Densitometer values 

Column 1-15: 

Card 3 (16F5.1) (DATAE) Densitometer values 

Repeat Card 1 and Card 2 and 3 formats for each successive group of 
3 data cards until ail desired data has been entered. End of data is 
indicated by a single card with some integer value other than 0 for 
Column 41. 


1 ) 
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PPOUdAM PBOMI lNPUTiOOTPUT.T*PF.'i«INPUT» UPC6«0UTPUT*PUNCM»r lUMPU 
Cjj.»a •?P 0<j«#H_t'.B(!P_I5.J3E5lei!»EU.j:0_0t»ELUH A DENSlTt-EUNCtlttN -OL-tHE 


'C ' " AZOLES !N THE iNPUt SCENP« BASED ON THE INTENSITY IN 4N 
C COPPESPONniNO DIPPPACriON PiITEPN neoative* the principal 
p C Pl'TPOr.lS.A.UISTBIDUllON OF ANGLES INFERRED, IN THE INPUT i - 

C IN lU RFOREE INCREMENTS FMOm 0 ID 90 UEOREESt 
C.t.t. INPUT 

_ £ IDEUll£XCAJ.lQN,JNtlUNi.JIDAf£l 

C TMHESHOLO VALUE IPERCTI 

Ip C BASF adjustment FOR APERTURE IHAOJI 

C MINIMUHJDIWEROFNCE TEST wAlUE (OTESTI 

(!i.^..outpu? 

c oriotnal data of intensities uff diffraction pattern 

C RASE value 

I A C ilAEtNUM deviation In' t'^'BaTa 

c avehade deviation 

£ D£VJAlIONS._ ^ . 

c ” diffraction pattern density function 

C SCtNE OEMSITY FUNCTION 

in D I mens I ON DAT AE I WI » PROPORIi£OJ-iUEJEiiai.AJ£MpiiaiiAA>l I SIASIiLTI Jig. 

c 


C...SET UP AXISES FOR NAPA ROUTINE 

ni't s'j»i*Te 

EA XAX|SU>«J«10. 

5_OlMLlNU£ — 

HUIOMPMOPORTION 
VTalOHUFAF ANGLE 

. . LTITm-IOMDENSIIr.EU 

3r LTTI IJlaJOHMCTION FO 

LlIIOIalOHR FOURIER 

LTt TUlatOHTECHM iOUE 

DO 6 JaStH 
6 LTlTUlalnH 

JS C 

C.iiERTFO MAJOR 'LTiDP ; " ‘ 

C 

Of 200 K*IiSO 

^ 

AP C...HFAD OAIA 

PFAD isiSiorTimjN.NUATEipERCT.BAOJ.OfEsf.MfH^ 

SIQ F')aHAriIA.I6«3F|0.S«ln 

IFIMIRAIL. OT. O) (J0_T0 300 - 

AS READ IStSOSI <DATAEtn«lalt|9l 

S»S rriRHAf MbFS.II 

c , •. . . . . 

C...i)ETFRtilNE BASE VALUF <LARrtPST DATA POIMTI 

r 

S'* OFi/IOIjO* 

HAAlaiiATAFU) 

Oit 10 .l«2«19 

IFCRASF.LT.OATAEUI ) RASPal)AlAEtJ> 
iO CitMlINIIF • 

5S C 

C...IFST FUR INCOMPATABIU IV BFTaEEN 3 AND lUO DEGREE REAUImOS. IF 

— - c Test is posttive'Then the data may be in error. 
c 

FlFSTlaABSIOATAFM >-«ATAP(lut ) 
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,.U!L 




13ri 


ia*; 


l4»n 




In 




PRlCCtlUlNO PAGE 


Mb n'»MAf(///,* ntVlATIONS«> 

-HHUe lA.l£SiJl_{J)E.V.(U J I»i 1 IS). . 

*(PITE(Atb3fc> 

AJO PnoMAl (///t* 

V(PjTEJ»>fWU) 

HPTTE(A»t»?0) 

670 P(IOHA) (///•• 


nlH'BACTlOw PATTbHN OENSiriT FU^ICTION VALUES*) 
ITEHPtl».I»H 19 ) 


DtNSITV FUNCTION VALUES FOP INPUT*) 

tHROeO(iT.aJ.I»a^JBL 

PUNCH TOb.f’RUNtNOATE 

PUNCH 7JO. (PROPORT) I) «t >1<10) 

Nl<MPT>|p _ . _ 

rt-A»»o. ~ ■ 

oil 90 J*i . IR 

( J ) . QT . YMA X ; YMAXaPR(Jt>OHl.tJ,) 

MA<«(rMA*».U*l0. 

VH#X»MAX 

VFA«»rH<iX/jO, 

call liHAPHIxAlIR. proport, YMAA.NUMPT.NRON fNOATE) 

T»5 rooiAMr*,I6) 

•TTTTTng HArtAMq.Al 

on TO 200 


?S0 WRITE (6,6|^«i) nrun 
6?b FT'RmST i*!*;* ' run‘ 
1 TO DEVELOP A DENSITY 

?.!!!) CONTI HUE 

300 CONTINUE 
STOP 

. ekd 


*,lR?**OiifcS NOT contain SlONIFiCANT OEVIATION 
FONCTTON*! 


SURROUT I NE graph ( X , Y , VMA a • Ni (HP r , NRUN , NO ATE ) 

DIHEN5ION XISOt.YlSq) 

LTlT)»10HPON NUMBER ' ~ 

LTiT2»5HDATE 

LABA-SHANGLE 

LAHY»1ohP«OPORTION 

call SET<.1,.B,.I,.8,0. , |80, ,0. <YMAX, 1 1 

. C*U.J?«(‘tTjSjA00»LAOY,J.0j4ill 

call PWRlf* 30 . 5 ,LAbX,S. 1 . 0 ) 
call PWRT» 382 ,HilS,LTm, 10 , 1 , 0 ) 

... . call HyM8P(NRUN,ZHj3l 

call PWRT I3'J5,'»60»LTlT2,R,i.O) 
call NURbflCflOATE,?rtI6> 

CALL D HDFmtI SHF 1 0 . 0 ,bHF 1 0 . 3 ) 

call HERlML(i8Vl7?«,ir 
call curve «x,y,numpt» 

CALL frame 
PF TURN 
EFO 




blank NOfHU*H| 
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Exhibit C-2. PROGflAM THETA. 

Program Name: THETA 
Subroutines Required: None 

Narrative: 

THETA is entirely self-contained. It derives a matrix of the 
three-space angles inferred by any pair of orthogonally projected 
angles from 2,5 to 87.5 degrees. This program outputs a punched 
deck of the matrix, which is utilized in program CONVOL. 

Control Card Input: None 
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t- 






in 


■1 I 




3^ 


«ii 


-*r 




S5 


C* • • • 

c 

c 

-C 


peoOflAM lH^T4CIMPUT«0UT9Ur<TAPt3BlNHUT<TAPE6«0UTBUT«PUNC*O 
.pAonHAH TwfTA IS Dt■SlQNe^ Tn oeTEHMlNfe THE THHEr SPACE AMPLE PnBMED Bv 
ALL PA[RS nF_PLANEg AnQlES IN FIVE UEONEE INCREMENTS SETfEEN 1 AND H9 

nlriBEts.fOP evamplE If a lInF fomms a in oEnhee anOLE In tme kf vie" and 
ihSM?! S c*Wfe ‘ iagg palcjlatp tmp 


C respective three space AN0LFI9.A UE6HELSI • 



^ COERTM (internal ASSl OrxMgWT OF CONVERSION FACTOH FOS OEOPEES TO 


PAUtANSI 
CNTnFU 1 RAOIAnS TO OEQRfESI 

-A-lT££T^>fl;.REPRES&NTlN3 0.90 Db WEES IN S DEBpeE In TERvalsIN THE kZ 


VtEiJT 

r V (VECTOR for ye VtEF> 

OUTPUT 


TO lAKPAV containing the three space ANGLES OF ALL PAIRS OF ORTHOGONAL 
_0.b!Ejy.si0M ■ V art).,y( i s i . to 1 1 a . mt ♦ th 1 1 ■ . i a i 


Co.ARStSN CONVFDSION factors 


CPE0M-.0I7A533 

cPT 0 Eo«sr.?v 5 Te 


§...<ier ' p nvE 6eaAee ii^eiiet^eNrs rnA' the k ano v mws 


P«F.S 


OP 23 lalilR 


“■(P*3.T‘ 

?5 continue 

00 33 lii.lB 

X(T)>X(Ii*COCi}TR 
V(Tl»TtI»*CDEnTR 
35 CONTI tjUE 


C...CSLCltLAte the three space angle associated kith with all pairs 



DO SO laltlR 
00 S<> JPltlS 
7^.TAN(X<H> 


VPa(TAN(A(Jll )/( TANIVdl I) 
C1«|.*Y2**? 


TO ( 1 1 j» » I TH ( I • JI «CRTOEOI 
So CONTINUE 


C...UENEPAfE output 


on luo IP1*I0 

PUNCH At»« (10<Ii.|)«Ja|«}B) 

..fijs.roRHAnpMiitii 


HnTTEIG,61lt) ITD(t«v'l tjal’IS) 

Aio Fm»NAf(H <1RF5.1I 

IPO CONTINUE 


STOP 

EM> 
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Exhibit C-3. PROGRAM CONVOL. 


Program Name: CONVOL 

Subroutines Required: None 

Narrative: 

CONVOL uses the joint probabilitir irthogonal pairs of 
planner projected angles to convo'' ' .ogonal distributions 

of angles. The process involves cycling through a matrix of three- 
space angles inferred by pairs of orthogonally projected angles. 

Each element of the matrix is assigned to a three-space angle inter- 
val, which in turn, is assigned the value of the joint probability 
of the orthogonal pairs occurance. These joint probabilities in 
each class are added, with the final sum indicating the probability 
of that three-space class’s occurrance. 

Control Card Input: 


Card 1 
Column 1-10 

(110) 

Card 2-55 
Column 1-80 

(8F10.5) 

Card 56 
Column il-20 

(F10.5) 

Column 21-30 

{F10.5) 

Card 57 
Column 11-20 

(F10.5) 

Column 21-30 

(F10.5) 


(N) Number of orthogonal 

pairs 


(THETA3D) 3-space matrix; THETA 
output 


(BASEX) Base value in D.P. 
(AVDEVX)' Average deviation in D.P. 


(BASEY) Same; Y view 
(AVDEVY) Same; Y view 
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TT 


1^ 


2i 




30 


pnTiH!i^A>l CoMVOL ■( fiif'‘Jr.'bdfPUT*TiPES»llMPGT«T*Pe6*nUTPUT*PJNCHl 

puni^HAH tn^VOL IS UFSlCNfO TO MEuee r «<0 0»TmOOOn*L PL*SEq olSTniBUTIONS 

C OP AMlil.Fb IHTO tHFIP TBOE ThWEE SPACE 0I^5TPl6UT I ON OF ANGLES# THE 

r TJBoCEuUftt TKVOLWrTjETTBHrfifNG fHE JOjNT ■pftOBiHti.l'fT OF ANY PAIR OF 

C Arr.LLSIFX.- ll< OEfiflEES IN ThE *1 ANO 2S OEOREES IN THE YZ PLANES! AND 

ASSIUNINO ’ (E JOINT probability TO THE ThR_EE SPACE *N8il_ ASSOC I ATEO, FIJM. 
,AL PAih‘ Tex 1- To* AND 25 DEGREE PAI'b FORMS A 9.4 DEGREE 


C 

c 

'C“ 

c 

c 

TT 

c 

c 

c 

c 

c 

c 

T7 

C 

C 

•*r" 

c 

c 


THF OHTHOe 

angle in three SPACE!, several COMSINATIDNS OF PLANER PAIRS CAN FORM 

pbhbahilIties of all the orthogonal pairs ahich form the Three space 

ANOLt. program theta OETERHJNES_THE_THHE£_SPjACE_4NiJLE_tOIWEO„Br.#U 

'rAlT,S'15F"ANGrE"S''rN five bfOpEE INCREMENTS OETMEEN 1 ANO B9 DEGREES. THE 
OLTpUr FROh PPOORAM THETA (MATRIX OF THREE SPACE ANGLES! ACTS AS INPUT 

TP t his Prog ram. . 

. I*Cpu i' 


N (NIlMRFfl OF orthogonal PAIRS TO BE MERGED! 

THETASO (MAT RIX O F THREE SPACE ANGLES! 

■RiSEx.SAfiEy^tHTSE^Alor ?5 f diffraction “ATTEPN-USEO BY FOLLOWING 
PROGRAM! 

AVUEWX.AVDFVY (AVERAGE OPVIATlON IN THE DIFFRACTION PATTE RN- USED 
TTv" THE ■pOurSw ITiI5"PROTRjMr 


XZ.YZ (ARRAYS CONTAINING THE PROBABILITY DISTRIBUTIONS OF THE 
ORTHOGONAL SCENES! 

:nuTPUT - - 


RROfld-lfl! (THREE SPACE lEAF ANGLE DISTRIBUTION! 
RR0B( 10-22! (BASEX.BASEY.AvOEVXiAVOEVY) 


-UTHENTrTTJN- tn??7T7rtBT35V2T 

dimension THETA30(19.19I .XZ(19I «VZ(19! tPn0B(22! 


35 




C...Rt An PROGRAM PARAMETERS 
C 

READ (5»560I N 


Srli Formaf (HO! 

DO 10 I-1.19 

10 READ (S.SIQ! (THETA3DI1. Jl >j»lil9! 


45 


C...SET UP heading for printer output 
-c 


WBITE(6»filu! 
610 FPPMAT(*1 5 

n — B5 — nr" 

?OEvr»f///| 


10 

-75~ 


15 

"BIT' 


20 

"B5 • 


25 

■90" 


30 


35 40 45 50 55 

HllSCX — iVDEV)T-«SET " 


6 


SO 


C ... RF A P INP U T IxrTn V D- T T-glSTfftHIlTTOR-OFmwOLES 'AND- TMElR'RESPECTTVT'BiSr ' 

C AnP AVERAGE DEVIATION VALUES) 




DTT IDOTTiTTN 

READ (S.5?n) (B(U.I!.I"l<2( 
READ (5.515! (0(L1 t I ! . I>1 > lol 

100 CPN1INUF 

515 FOrMAT(8F|0.5) 

520 FORMAT (10X.2F10, 5) 

D O 20 P LI pTTN 


on 125 1*1.19 
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Exhibit C-4. 


PROGRAM ORTHOG, 



Program Name: ORTHOG 
Subroutines Required: PLOT 


Narratl ve: 

ORTHOG determines the three-space distribution of angles from 


pairs of digitized orthogonal photos. In addition, it outputs a 
microfilm plot of each view. 


Control Card Input; 


Card 1 
Column 1-40 

Card 2 
Column 1-10 
Column 11-20 
Column 21-30 


Card 3-n 
Column T- 10 

Column 11-20 
Column 21-30 


(4A10) (LTIT) 

Title for plot 

(110) 1 

NPLANT) 

Plant numbering sequence 

(110) 

.NLEAF) 

Leaf numbering sequence 

(110) 

;npts) 

Number of digitized 


points/ leaf 

(F10.5) (Z) 

Z coordinate of a digi- 

tized point 

(F10.5) 

[X) 

X coordinate 

(F10.5) 

(Y) 

Y coordinate 


continued for as 


many colimins and cards needed to describe a leaf 


Cards n-+ 
Repeat cards 


number 3 and 3-n format for remaining leaves and plants. 
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lu 


15 


ett 


2S 


^ir 


3 S 


♦0 


45 


‘il) 


55 


PfiOORAM OHTHOOt Input »OUTPIll • tapes- input tT*Pt6«0UTPUTt punch tFlLHPtl 

PROGRAM OHTMOO IS PESIGNCO TO TAKE THE DlGIT Ig EP OATA F ROM ORTH OGONAL 

— c vZj Al^'O COnsTPUST aN aBsTPaCT TnPCE SPAtf ^EPpEsEnTATION 

C Of THE ORIGINAL PLANT, THE PROORAH SAMPLES THE THREE SPACE LINEAR 

c segments comprising the plant to oev elqp • pistrib utIon o r angles 

lIaves of the plant. 

c input 

C NPLANT {BOOKKEEPING NUM BER O F THE P LANT I 

^ •"■■ULE*^“(BOOKKeEPIN6 NUMBER ASSIGNED TO EACH LEAF OF THE PLANT » 

C NPTS (NUMBER OF POINTS DEFINING EACH LEAF) 

C LTIt (ARRAY containing THE TITLE FOR TME_MICMniJt PLOT^I 

C ? (V£c1r£»H COfTfAlNlMG THl Z cDOROfNATES) 

C A I (VECTOR containing THE X COORDINATES) 

C X{ (VECTOR containing THE Y COUPDINATESI 

- -C Ll TlfrtERNAtVY'A^rPNEp vARttBLE INOlciTING THE NUMBER OF LEAVES) 

C MTEST (AN INTERNAL VARIABLE USED TO FLAO A NEW MICROFILM PLOT) 

C PLTX»PLTY (ARRAYS USED TO HOUSE T HE VALUES TO BE PLOTTED) 

-T - prod (ARRAY' 5f" THE PflOBSB I L f TIES FOR ANGLES 0-90 DEGREES) 

C COSTH (THREE SPACE ANGLE FORMED By A LINE SEGMENT) 

C OUTPUT^ 

C ■ PflOB '(blSTRIHUTION OF THREE SPACE ANGLES* PRINTER* PUNCH) 

C FILM PLOT OF THE XZ AND YZ VIEWS OF THE PLANT 

C 

■— "C 
C 

DIMENSIO N Z ( 5*IS ) ♦YIS*IS) *X(S*15) »PRgfljie) *KK(5) *LTIT (A) 

blHENSTON PLTX()8TVPLfY(lB) *PxZjl5)tPYZ(rffr 
C 

C,,. initia lize PARAM ETERS 

C 

H(EST-1 S NLlNESaO 9 Ll-0 S RSUH-0. S RXSUM.n. S RVSUM-0. 

00 3 I-ltlB 

PLTX ( irio 

PLTYtI)-0. 

PAZ(t)-0. 

P*Z(i)-0. 

3 pRop(n«o. 
c 

c;7;PEJir-pon-TTTi;r 

c 

READ (5*500) (LTIT(I)*I*J*A) 

SOO FORMAt(AAlO) 

C 

C...PEAU variables DESCRIBING A LEAF 

■ c ■ ■ 

10 READ (5*5201 NPLANT.NLEAF.NPTS 
520 FORMAT (31 10) 

C ■ ■ ■ ' ■■ ■" - 

C...TE5T FOR A NEW PLANT 
C 

11 IF(NPLANT.6T,MTtST) 0(1 TO 200 
C 

C...TEST FORENU OF DATA SET 

C 

)2 IF(NL£AF.EO.O) go TO 1000 
C 

■ 'X',,; ASSIGN internal PApAmFTERS 
L 

L1-L1*1 





147 



KKlUl'«NPtS 

COOROTN»T€S OF i LE*F ! 

C 

BUO (i*sasi (l(H*lli*(Ll*nfV<Lltll«l«l*NPTSI_ 

— 5?5 FOUKTriSFro".?} ’ ■ 

c 

C,,.C*LCULATt THE THREE SPACE ANOLES FOB ALI, LINE SEQHEN TS DESCB 1 B1N8 h L E>F 

C ■ ■ 

DU too L2«2.NPTS 

0X«)l|Ut«L2)-)l(Ll<L9>l) _ 

6TiiriLl«L2)-viU •L2-l> 

OZ>ZtLl'LE)-Z(LItl2>tt 

B»SORT(OX*»2*Oy»«?*OZ*»2| 

* ' UZ I /B 

CUSTHb 90.« t ACOS (COSTHt •Sr.ZBSrTBOlSl } 

RAZ»S0BT 10X» »2*0Z»«Z) 

THiiZaA&§|OZr/RFZ 

ThXZ»AC0SITM*Z»»3T.Z95TT95l3l 

THXZa40«-THXZ 

Rvz;sflffnef«52?»l««z> 

THYZ«*6S(0Zt/RVZ 

THTZaACOS t TMVZt *57,2957705131 

TH?ZiWi^?Z 

C 

C.. .ASSIG N EA CH THREE SPACE A NGLE TO OWE OF THE 5 DEOREE IN TERVALS BETWEEN 

“C ■# XN6 90 DIQSCER; fhg LENGTH oF T^E LINE SEOHENT INDICATES THE RELATIVE 
C vetOHTiNO OF THE ANOLE 



toll 


65 


TO 


76 


1 , 0 “ “ 


ttS 
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1<?U 


1»S 


l3U 


1 J5 


1>»U 


U5 


ISU 


139 


l6U 


145 


17U 


1/5 


ifUO CUNTlNUe 

MieST*H>'l«NT 
SOMmO. 

XSUMbO* 

YSUhcQ. 

DU aio Ui.iB 
X5UH«XSUH»PX2m 
YSUMiiYSUM*PVZ<II 

aid 5UMBSuM*PRun(ii 
DU ?zo i=j»ie 
PXZ(I)»P*ZtI|/X5UP 
PYZ(I)«Pif2Hl/YSU«> 
azo pDuRdi «ppun<n/suM 
c 

' C..,SFNt«ftT PHiNTER OUTPi.lf 
C 

NP«NPLANTrl 

WPITei5«62S| NP 

(.as F(iRH*n///i* PMne*fliuTY*t/»* of AN6te»t3Xi« s in i? ao ? 

IS 3d _ 35. 4U 45 so 55 _00. 65 70_ 76 50- 55 SO 

a FOR PLAHT*tl3«/| 

NRnEl6«(>30l <PPOR(I) ilaltlH) 

630 furraxi* ^a^s^Ace •ijLBE.s.ej 

WRITE I6«6S51 tPX2<li.I»t»rat 
655 FORMAT!* XZ ••18F5.21 

WH1TE(6 «6651 lP YZMl»l«l.lttl 

665 FORMAT (• YZ ••iflFS.21 

WHlTE!6t6351 RSUM 

63S_ FURHAT 1 // 1 * t o tal P L*HT L eNO TMi.». F6.3> 

C 

c...0ENbRATC Punched output 

c 

BASE* I . 

PUNCH 705»NP.flSUM 

F05 FURMAI(H.O.»F10.5l 

PUNCH MO» (PROnm .T«1«1A) 

710 FORMaT(8F1u,S» 

0 PUNCH 707tRX 5UM 

707 format ( lox.aFin.s) 

C PUNCH 710. (PXZII) tlsi .161 

C PUNCH TOIiRYSUM^- . _ 

C PUNCH 710. (PVZ(l) ,1«1,)6) 

C 

C... generate FILMP1.0T OUTPUT 
C 

DU 165 Isl.Ll 
«T£ST»I 
Mi>KK!il 
DO 160 J-l.HI 
PLTX(J1«X(I.J) 

PLTY!J1»ZII.JI 

160 continue 

C CAUL PuOTCPLTX.PLTY.I.KTESr.Ml.UTITl 
165 CUNMNUE 

(lU 175 1>1.L1 
KIESTsl 
Hl>KK(ll 
DO 170 J«1.M1 
PLTX(J)=Y1I.JI 
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ISO — 


las 


l7J>_CONTlNUt 

C "cAtL ■^5Tl>LTK.MLTYt2,l(TESTtMirLTIT> 
ITS CONTINUE 
C 

C'iiVBElNrri»Llze PARANETERS 

c 


A"°J*|-jNgS« 0 $ n«0 % RSUH« 0.1 RXSUHaO. t BYSUM.B- 
DO Z50 I«if 18 * 

PYZ(I)«Q> 


PTZm«o 


14B 


‘ asiTPHOa 


(in»or 


00 TO 12 
__IOgO CONTINUE 
STOP 
END 


to 


subroutine PLOT(X,Y«MVieWtKTESTt«Kfl,TlTl 
DIMENSI ON XT ZOt ♦Y|2«t .LTITTBI 

^F^(fTEST;I^T7Ir ‘TO t(j'*s " - - -- 

call frame 

CALL 5ETTO«tl»tO.»|»«.».,».«»i,,7.,i| 
C*Ll StTLil^E(i) 

CALL FRSTpTt,3.5,..5t 
call STMINOtLTlT l 
If (NVlEU.EQ.2l Qo TO 20 


IS 


20 


2S 


lViem.iomiux views. 

so TO 30 

20 LVIEN>1uhsuy Vitus. 

30 call FR5TPT|3.,-.Sj 

ca ll STHlNQtLVIE WI 

CAlL LINE7>3<«0'.>*1.i6«I 
call LiNe(-.i.,*,.3.,ft., 

CALL L INE ( J. 1 6. 1 3 . ♦ 0 . I 

call LTfiE73.,0.t-'3..0.» 

40 M2*KK*1 
00 100 lal.HZ 

xPTi.xTTi 

YPTl.Ydt 

XPT2»XI1»1| 

yPt2^y(i»h ' • — -- 

CALL LINElXPTl*YPT|tXPT2tVPT2) 
100 CONTINUE 

RETURN '■ 

END 


^ 4 . 
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Exhibit C-5. PROGRAM SUFU. 


Program Name: SUFU 

Subroutines Required: MODIFY ADDSUB 

FEATURE TRNSPO 

SMEANS MATMPY 

INVERT 


Narrative: 

SUFU is a classification routine for comparing a spectral signa- 
ture of an area with unknov'^n LAI with the spectral signatures of 
several areas of known LAI, in order to infer the LAI of the unknown 
area. Subroutines MODIFY, SMEANS and FEATURE are user defined and 
allow for the modification of the input data. Output consists of a 
C-factor matrix, calculated divergence, and a printer plot of diver- 
gence. 

Control Card Input: 


uara x 
Column 

1-10 

(110) 

(NSIG) 

No. of "effects" 

Column 

11-20 

(110) 

(NPLT) 

No. of "plots" 

Col umn 

21-30 

(110) 

(NWAVE) 

No. of v/avelengths 

Column 

31-40 

(110) 

(MODVAR) 

Flag to modify co- 
variance matrix 

Column 41-50 

(110) 

(MO DM) 

Flag to modify mean 
vector 

Col umn 

51-60 

(110) 

(IFEATV) 

Flag to derive fea- 
ture vectors 

Card 2 
Col umn 

1 

O 

(AlO) 

(I DENT) 

Identification of 
"plot"/ "effect" 

Card 3 
Column 

Card 4- 

1-10 

-n 

(F10.5) 

(MEANS) 

Mean vector of "plot"/ 
"effect" 

Col umn 

1-80 

(8F10.5) 

(CO VAR) 

Covariance matrix of 
"plot'V'effect" 
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Exhibit C>5. Continued 
Cards n-+ 

Repeat cards number 2, through 4-n for remaining mean vectors and 
covariance matrices of all of the "effects" and "plots". The order 
of entry is all of the "effects" data first* followed by all of the 
"plot" data. 
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5 


10 


IS 


20 


as 


30 


35 


4e 




SO 


55 


PHOGRAM SUFU (INPUT. OUTPUTiTAPE5>lNPUT»TAPC6itOUTPUTi 

Ol HENS iO N B(A.41 .CMAT<* «A> 

TIHE^TOKI CUVAH<4t4il3|.I0£NTTl3i *SIGl (4,4| ,St02(4.4| 

DIMENSION XAVKA.l) •XAV2(4.11 .CUtA) tV(4.i| .01 1 .4} «Z (4. 1 1 

DI MENSION CFACTl O .l OltCPACTa 110.31 . TEMPI*. 41 .SUrUDd O) 

OlHENStliN Wl(4.4)iHa'|4.4|.XAXISU5» 

DIMENSION EI4.II«P(4.U.FACl(4i4».PACa(4.*i 

HtAL MtAHS(1 3t41 

C 

C.,,HEAU PROGHAM parameters AND DATA 

c 

READ (5.505) NSIO. NPLT* NMAVE. MOOVAP. HUDM. IFEATV 
5u5 format (GUO) 

KSTPaNS 10 ♦ NPLT 

BiTlOU (.alfKSTP 
READ <5.5101 IOENTIU 

510 FORMAT (A1 0> 

READ <&»iil5l (HEANS(E.J) .J»1.NUAVE> 

OU aOU 1*1 .NMAVE 

READ (5.S1SI (C0VAR(I.0.L1»J»1.N MAVE> 

515 FURHA1 ibl'lO.Si 
aOO CUNTINUE 


irrj.PRINT INRUT 

t 

MRITE (6.600) 

“ ' ■" MHnE (6.001 )'NS1 o 7 NPUfTSvAVE. MODVAR. MODM. IFEATV 
eol FURMAT<* NSI0*4.14./.* NPLT**. 14./.* NNAVEa*.I4./. 

I* M00VAR**.I4./.* M0DH**.I4./.* IFEATV«*.l4. / /// /) 

Dims 1«1.KSTP 
HRITE(6.63a) IDENT(l) 

632 FORMAT I /////.* RUN — — 

HNITE (6.6^7) (MEANS (l.j) .0*1. NMAVE) 

637 format (//.36X.* CHI CH2 LH3 CM* CHS 

1 CN6 CH7 C HO*.//.» MEANS B Y CHANNEL*. 15X.BF1 0.5 ) 

HH1TC(6(647) “ 

64T FURHAI|/«'.* COVARIANCE MATRIX*) 

00 255 X* I .NMAVE 

255 MR1TCI6.642) (COVAR(K.J.I) .J*1.NMAVE) 

642 FORMAT liSXtOF 10.5) 

215 CUNTINUE 

BRnnoTWoT 

KST*NS10 ♦ 1 

IMMOUH.eo.l) 60 TO 216 

ooTO'sr? 

c 

C... SMOOTH MEANS 



216 CONTINUE 

CALL SHEANS (NSIO. NWAVE. MEANS) 

“iRTTETStiSUn 

618 F0RMATt*l*.5X.* SMOOTHED MEANS USEO IN CALCULAT IONS ....*»// I 
Oil 21H laltNSIG 

— smffiTrrs.ii?) i,iHCARs7i;;rr;j=nww9n ~ 

617 FURHAT(lUXtI2tSX.0F>0.5) 

217 IF (IFEATV. EQ.l) 80 TO 7«3 _ 

00 T0"^19 

C 

C.. .DETERMINE NEW FEATURE VECTOHS 


oRIdlNAI, PAGE IB 
jF poor QUAUTU 
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bO 


6S 


TO 


lb 


■W 


as 


90 




loo 


105 


"nr 


us 


Tq3 call FtATUBElNUM tWE4NS.MUL ltHatLa.M3.LJ , H« « L« I 
WHITEI6*eo2» NUM* H)iLltM2tL2i)l9«t3fM4tL« 

802 F0HM*Tt/////»I2t* NEV FEATURE VECTORS Otfit OEVELOREO. THEV ARE.,.. 

BAWfl2«o / BAND0«12M 


c... enter plot loop 

C 


219 DO 400 LaKSTiKSIP 


_itA.l CENTlFY appropriat e PEAN VECTOR ANO COWAwIANCE WATRIK FOR PLOT 
DU 120 t-l»NHAVE 

DU 220 JaliNNAVE 


TEHPIltJIaO. 

220 CUNTINUL 

IF (WOOVAR.EQ.il 80 TO 230 


DU 225 laltNNAVT" 
XAVlIltUaNEANSILtU 
DU 225 JattNMAVE 
TbMPUtJUieOVAAil.J.Ll 


225 CUNTlNUb 
OU TO 2*5 


C... MODIFY WEAN VECTOR ANO/OR COVARIANCE MATRIX 
C 

-"“niTCDBTnraE 


DO 233 laltNWAVE 
232 XAVtllUlaMEANSlLiI) 

DU 231 laUNbAVE 

00 231 JbIiNHAVE 

231 CM AT 1 1 • J|*COVAR(It J t LJ 

OU 240 latiNNAVE* 

DO 240 JaltNWAVE 

fCTPUiJ|BB(Itj} ~ ' 

240 CONTINUE 

VHITE(C>H>8H IDENT IL) 

68 rFoRMirW/,* 5fo6 1 fTEo mEan VECTOR ano/or covariance matrix foF 
1 MUN* » ■ AlO) 


C...PfllN| MEAN VECTOR ANO COVArIaNcE MaYRIx OF PLOT USED tW CALCULATIONS' 
245 VRITE(6f715) IXAvl (J«I) . Jal,*| 

ns FUHHATIUPlo.S) ' 


DU 10 lal.A 
MR1TE(6.715> (TEMP(1,JI,U*U4| 
iO CUNTINUE 


C... ENTER EFFECT LOOP 

DU 400 Hal.MSlO 

c 

^...IUCTUfy appropriate mean VECTOH *NireOvARt»NCE HAT ff I k POWEFrECT' 

C 

DU 247 lal.NHAVE 


5iT 24T Ja >, .shave 
51Cc(UJI*0. 

247 CUNTINUL 
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_UQ_ 

IT IMOOVAR.EO.I) go to 260 
DU gqii l»l.wy*up 


XAV2(lit)«HEANS<HiI| 
00 2S0 <^BltNHAVE 

125 

250 CONTINUt 
00 TO 2>0 
C 



c .. .modify mean vector and/or covariance MATRIX 



260 CONTINUE 


130 

00 262 iBltNNAVE 
262 XAvEIlf DsMEANSlMtll 
DO 261 lal.NNAVE 


135 

DO 261 JaltNNAVE ' — 

26) CMAT(l«J|aCOVARItiJtM| 

CALL MOUIFYICMAT.B) 


00 205 iBliNHAVE 

DO 265 jBltNUAVE 
S102llt JtsBtltJ) 

14U 

265 CONTINUE 

NHITEI6>601| IDENT(M) 
C 



C... PRINT mean VECTOR AND COVARIANCE MATRIX OF EFFECT" 

USED IN CALCULATIONS 

...... 

2T0 HKITE<6t7I5l IXAV2(j,i) i.)b1iAi 


145 

00 20 I»lt4 

HRITE(6iT15) (SI02(1«J|»J81,4) 
20 CUNTINUE 


_iS9 

c 

C... CALCULATIONS FOR C.FACTORS 
C 



NRITE(6i6S6) L*M 

656 FORMATI///*# C-FACTORS FOR PL0T»iI3»« AND LAI 

ialculation steps completed. ..•) 

EFFECT4.13,/.* c 

155 

c... DETERMINE INVERSE OF COVARIANCE MATRICES 




WPITE<6f651) INHERE 

651 format l/yillO) 


C., .assign COVARIANCE MATRIX OF PLOT TO NORKlNG MATRIX 


165 

DU 2 7 1 a 1 f 4 
00 27 Jal,4 
S101|ltJ).T£MP<I,J| 


LUNTlNUt — 

t call invert ISiei.NWAVEl 

call HAtftlX(10t4,*tOt5IGlt»iOETRHl 

iro 

JUUV 

tOUO HHITE(6*69n (SlOld.JI.Jal.A) 
691 FORMAT I4F20.S) 


1KMEREB2 

NHlTE(6t6Sll inhere 

C call invert ISI02,NNAvE) 

CALL — 

OU 1001 Ial«4 

1001 NRlTE(6t691l (Sl62 1 1 . J) f J.l .4 1 
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•( ■ j 


ISO 

c 

IWHERE*3 

HKITE«bi65n inhere 


las 

CALL AOUSUB |AAVIt!;AV2fE«NNAVEtItVl 
NHITE(6t692» IV(Iil»*I*li«l 

602 FORMAT IF20.5«3I/.F20. 511 


inhere** 

WHITE (btbSl) IWHERE 

call Auusue ix»vatK*tfi»r,Nw*vEjl«zi 

190 

HHITE(Ot692> (Zllil>*I*l«41 

c 

C...TAKL TRANSPOSE OF SUBTRACTED KEAN VECTORS 


c 

IwHERE>9 

HRlTE46t6Sll INHERE 

195 

CALL TRNSPO |r<CtNHAvE«ll 
HRITEI6«691I (C(l<JltJ*li4) 
INHERE** 


200 

nh1te(6i65iV inhere 

CALL TRHSPO (ZtOf NWAvEtl 1 
HHITEI*f691) (D(l«J)»J>lt*l 




C...HULT1PY UERIVEO VECTORS 
C 

205 

IWHERE67 

NKITEr6f6Sl> INHERE 
CALL HATHPV (T»C«Hlt*tll4l 


OU 1002 l*ltA 

1002 MRlTE<6i69t) (Hi (I i Jl t J*1 •*! 
lWHERfc.»B 

210 

HHITE(6*6S1I INHERE 
call HATHPV |Z«0fN2«4.1i4| 
DU 1003I*1<4 



lOOi NRlTEIbi691) (N2 1 1 f 0) • J*1 *41 


215 

C...HULTIPLY INVERSE OF COVARIANCE MATRICES By (HE DERIVED MATRICES 


c 

IWHERE*9 

WHITE l6i6Sl) INHERE 

220 

CALL HATHPV |SIOlfHl«FKCTiAiA>A) 
DU I004'I*I«4 

1004 NHlTE(6i691l (FACl ( I t J) i J*1 |4> 


IwHERE-10 

UKITE(6<6Sl> INHERE 

CALL HATHPV (SI02tH2«PAC2i*t4*4» 

225 

UO^IOOB 1*1 14 

1005 HRITE(6t691) (FAC2 ( I f J> » J*1 •*) 
C 



C*..TAKb THE tRACE or THE FIRAL CACCOCATEO mATrICes 
C 

CUlaQ. 


230 


C02>n. 

DU 325 J*1«NHAVE 
cuts CUl * FAClIJtJI 

235 

C02* CT2* TACTI J7JT 
325 CONTINUE 
IwHERE-il 
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d 40 


24 !> 


asD 


25S 


SbU 


Wt'ITEtbibSJl IWHE‘»E 
LS>MFT=L - NSIS _ _ 
CCAcfl (LSHF1,M)«C0l 
Cf ACT2IM»LSHFT)»C0a 


666 FOR^t ^C^ACTOh’fM^ I-*I J!,!!!/ w? 

1 I 4x.>iri>nM Lai FFrECT»>13<* TO PLOT*.i3** IS^tFgO .5 1 ////////'' 1 


400 CONTINUt ._ 

L 

C,, .PRINT C-FACTOR MATRICES 


"■ WHITE (bf 690) 

WH1TE<6>6S2) 

652 FORMAT (/////.* j:-FACTflHJ|}MEl£feS,ijLa?j/Fi 
WHITEI6*6B2) 

662 FOHMAH///** FROM PLOTS TO EFFECTS*) 

00 420 laltNPLT __ 

W«ITE(6»"671) iCFACTl (1*J) |J"1|N5I0) 

671 FORMAT(10F12,3I 

j420 CONT_lN>^ 

WHiTE»6t683l 

683 FORMAT (///«• FROM EFFECTS TO PLOTS*) 

OO 425 l»i,NSIO 

WHI TE (6f67U (CFACf2U*J) tJ*l*NPLT) 


425 CONTINUE 

00 601 ialiNSlO 
601 XAXIsrtlal 

00 600 Lsl.NPLT 


p — - - . . ..- ... 

^',T:ALtULATE SU-FU DTsTaNC"ES FROM" (>F ACTOR HATMICES 

t 

WWITei6t69Q) _ - - 

690“fORMATI*1*) 

270 OlSTslUOOOOO. 

mN«l ■ 

DO 500M=1.NSI6 

SUFU0(M)«(S0RT(CFACT2(M,L)*CFACTI IL.M) ) |/<SQRT (CFACT2(M,L) ) 

I*SORT tCFAai.(L»M) ) ) 

275 ' 17 (SUFUOIM) .LE.OIST) GO TO 40 

GO TO 60 

40 D1ST=SUFuU(M) . .. 

CTfliM 

50 CONTINUE 

2BU 0 _ ^ ... - 

■ Ci.. OUTPUT RESULTS 
C 

WHITEI6i 662) L.M.SUFUDIH) . . rj.- , . 

66^ FTfflKm* T>ir^£P*HA0rLlT?‘B£rWEEfii PLUT*.l3t* and SlMULTfED LAI 

285 IEFFECT**1J.» IS*.F10,5I 

500 continue 

■ - WHlTET5T6?7rtiISf.ieENTtMlNV 

667 format I///.* THE MINIMUM SWAIN-FO CLUSTER SEPARABILITY DISTANCE 
1 IS».FI0,5»* ASSOCIATED WITH LAI EFFECT'.AtOI 

-2no CALL RaPIT iS'.XAXIS.SUFUO.". fNSiG.HL.HH.VL.VH.IOrt EFFECT . 

HUM OITERG. .lOHOtVG. PI.UTH) 

600 CONTINUE _ . .. . 

- ^ - 

END 
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s 


10 




20 


subroutine feature (NUMtXAVOt 

OIMENSIUN XAVOaStAliXWORKllStA) 

» L1«2 


HliLUH2iL2tM3iL3fHAiUA) 


Mt'ca s U2s0 

H3»L a U3»0 
M <*=0 i LArO 


nu 10 J>1(4 
DU 10 1«1«I3 
XI» 0 RK( 1 >J>« 1 , 

10 CUNTINUE 

DO 20 lilti a 

XWORK 1 1 1 lla' XAV6 ( t tMl I 
20 CUNTINUt 
NUHa 1 

DU 30 JaltA 
DU 30 i«lfl3 


XAVOaiJ»«XWORKCI.J) 

3crcDffritJ0E 


/ 


A«VO<ItLll 


RETURN 

END 


5 


subroutine SMEANS (NSIGiNIiAVE (MEANS} 
REAL MEANS (13(A) 

ISTOPr NSlO - I ' 

DU 217 Jal(NUAVE 
DU 217 i.i.lsTOP 


2 1 7 CoS nrii^ n-1 . J) /A. ) ♦ ( MEANS ( I ( J1 /2 . | (ImW NM H 1 , J ) / A , I 
_RETUflN 

END 


SUBROUTINE modify (C,B) 

DI MEN S ION 6 (A (A>. C(A .AI 

OU 10 1«1(A 
DU 10 Jal(A 



IF(i.EQ«Jl Bir(J}a,0000r 
10 CONTINUE 
REJURM_ 

END ' ■■ ■■■ 


S3RIGINAI] PAGE IS 
OP POOR QUALITY 
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SUUHOUTINE INVERT (A,N) 


c , 

5 

c AuBWOuilNE INVERT USES THt OAUSS“J0R0AN EUIHINATION METHOD TO 

C HtPLAtt THE MATRIX. A. 5ITH ITS INVERSE. 1 ME DIMENSION STATEMENT IN THE 

C MAIN MHUGHAM MUST HAVE A LISTED MATRIXt fcUXbOZ HQtlEVERj._n<6 ACTUAL -NUMUEB- 


C MAIN PRU6HAM MUST HAVE • MATRIX UlMtNSlUNbU mHICH IS tUUAL 
C than THt 50HKIN6 MATRIX IN THIS PROORAM. I HE ACTUAL NUM8ER 
C AND COLUMNS Of .A. THAT ARE U^ED IS SPELlflED fl¥ N. 

OH 

OF ROWS 

10 

C 

DIHENSIUN A|4tA),S(A)«CI4)tLZ(A) 
DU ID JBl.N 


ID LZIJIbJ 

DU 20 iBliN 

.. 

15 

TbA( 1.11 
L»I*1 

LP»I»1 



If (N-LP) 14.11.11 
11 DU 13 J-LP.N 
VBAII.JI 


■ 1MABS(V>-ABS(V1 1 13«13«12 

12 K«J 

25 

13 CONTINUE 

W«IT£(6.5151 Y 

AIS FURMAT (HF10.5I _ 


1 * 6u is ufifN 

CU)>A(J«K) 

A(J>K>bA(J«I) 

3fl 

A(J.I1«-C1J)/Y 
A(I«J)sA|I.J|/Y 
15 HU)*All.U) 


35 

'A(l>licl.p/Y 

J=LZ(I) 

L2(1)bL2(K) 


l/Tk7T3 
DU 19 K»1|N 

IE<1-KI I6il9«16 ... - 

«0 

16 DU la JBl.N 
IMl-Jl 17.10.1? 

17 A(K.J)bAIK.J>-B!J)*C|K) 


i'4 CuNTINUt 
19 CUNTINUt 
2D CUNTINUt 

-*5 

DU 2UU lol.N 
If(I-LZtll) 100.200.100 
loo K*l*l 


50 

DO 50(1 JbK.N 
Ifd'LZUn 500.600. 500 
6U0 MbLZ(I) 


■ LZ{M«LZijy ■ ■ 

L2IJ|bM 
DU 7UV L>ltN 

55 

— -tfu.Air.-u 

aiI.LIbA(J.L) 
7U0 A(J.L><°0<LI 



2U0 CUNTINUt 
RtTURN 
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subroutine AOUSUB IA.B«CiM«NiYI 


c 

5 

0 . 

c 

u 

>I...SUC 9 R 0 U 11 NE AOUSUB STORES THE RESULT OP 
result UP A-B IN V, MATRICES AiBtCtANO 
MAN matrices in THE CALLINO PROORAM. 

A*b IN C and THE 
r HUST ALL BE DIMENSIONED AS 


c 

dimension a IH tNI fB(H|N> iC (HtN| «V(M,N» 
DU 1 IsltH 


lu 


OU 1 JcliN 

CIItJ)sA(l,J)*BII<J 3 
1 YIJtJ|cA(I,J)-B(!tJI 



RETURN 

ENU 


subroutine MATRPY (AiBtC.MiNfLl 

c 


t.,,..SUBROUTlNE MATHPY MULTIPLIES TWO HATRICtS <A«Bt TOGETHER 
C AND STOKES THE RESULT IN C. THE UlMENSluN STATEMENT IN THE 

C main PHOORAM MUST HA y E A OIHENSIONE Q AS AN MAN MATRIX. 

C B MUST bE OIHC.NSIONED NXL ANU C MUST MXL# 

C 

_DiMENS10N tBiNtLI «C(M»U 

00 102 lelfM 
10 OU 102 Jal.L 

S» 0 . 

UO 101 KalfN ■ 


lOI S>S*AI lfKI«BtKiJ) 

1 0 2 C<liJ)«S 

1 B RETuBn 

CNO 



c 

subroutine TRNSPQ (A,BtH,N) 


S 

c* 

c 

c 

....THIS SUBROUTINE CALCULATES THE 
AND STORES THE RESULT IN b. A 
MAIN PHUORAH AND MUST BE NAM. 

TRANSPOSt UP THE MAN MATRIX A 
MUST BE UlHENSIONEO MAN IN THE 


c 

DIMENSION AlM.NI «B(N.MI 
OU 10 t«l.M 




OU 10 Usl,N 

' 

10 


10 BIJ.UMAd.j) 
RETURN 


END ■ ■ • 


■ -fcS. S 


IS 



